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ABSTRACT 

We study the uniqueness and robustness of acoustic signatures in the cosmic microwave 
background by allowing for the possibility that they are generated by some as yet unknown 
source of gravitational perturbations. The acoustic pattern of peak locations and relative heights 
predicted by the standard inflationary cold dark matter model is essentially unique and its 
confirmation would have deep implications for the causal structure of the early universe. A 
generic pattern for isocurvaturc initial conditions arises due to backreaction effects but is not 
robust to exotic source behavior inside the horizon. If present, the acoustic pattern contains 
unambiguous information on the curvature of the universe even in the general case. By classifying 
the behavior of the unknown source, we determine the minimal observations necessary for robust 
constraints on the curvature. The diffusion damping scale provides an entirely model independent 
cornerstone upon which to build such a measurement. The peak spacing, if regular, supplies a 
precision test. 
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1. Introduction 

Cosmic microwave background (CMB) anisotropies provide a unique window into the early universe 
through which we receive information on both the model for structure formation in the universe 
(e.g. Efstathiou, Bond & White 1992, Gorski et al. 1995, Ostrikcr & Stcinhardt 1995, Crittenden & 
Turok 1995, Albrecht et al. 1995, Durrer et al. 1996, Dodelson, Gates & Stebbins 1995) and the 
background cosmology (e.g. Bond et al. 1994, Hu, Sugiyama & Silk 1995, Jungman et al. 1995, Seljak 
1994, Scott & White 1995). In the simplest models for structure formation, based on the gravitational 
instability of initial density perturbations, the acoustic signature in the anisotropy power spectrum provides 
a clean and unambiguous means of measuring all of the parameters of a Friedman-Robertson- Walker 
cosmology. In particular, it offers a standard ruler with which to make a classical test for curvature in 
the universe (Doroshkcvich, Zcl'dovich & Sunyaev 1978, Sugiyama & Gouda 1992, Kamionkowski, Spcrgcl 
& Sugiyama 1994, Hu & White 1996). However, structure formation may have proceeded by a more 
complicated route. Recent investigations have begun to probe the acoustic signature in texture (Crittenden 
& Turok 1995, Durrer et al. 1996) and string (Albrecht et al. 1995, Magueijo et al. 1996) models for structure 
formation. Their predictions differ strongly from the standard inflationary case and suggest that perhaps 
without prior knowledge of the correct model for structure formation, the information contained in the 
acoustic signature cannot be extracted. In this paper, we focus on two general questions: does the acoustic 
signature uniquely specify the model for structure formation? how robust is a measurement of the curvature 
to changes in the underlying model? 

As discussed in the appendix, the question of uniqueness is especially interesting in the case of the 
standard inflationary paradigm. Inflation is the unique causal mechanism for generating correlated curvature 
perturbations above the horizon (Liddle 1995). All other causal mechanisms generate significant curvature 
perturbations only near horizon crossing. We will refer to these alternate possibilities as isocurvature models. 
A unique signature of superhorizon curvature fluctuations can be used as a test of inflation. By generalizing 
the external source formalism of Hu & Sugiyama (1995a, b; hereafter HSa, HSb) to include backreaction 
effects in §2, we find that a large class of isocurvature models carry a distinct acoustic signature that can 
be easily distinguished from the inflationary case, independent of the curvature and other cosmological 
parameters. The distinction lies in the gross properties of the spectrum, not small or subtle shifts in the 
peaks heights which require very high resolution measurements to see. As shown in §3, tell-tale features 
include the harmonic series of peak locations, the alternating relative peak heights due to baryon drag, and 
the diffusion damping tail. From this study, we conclude in §4 that the standard inflationary model with big 
bang nucleosynthesis (BBN) baryon-to-photon ratio bears an essentially unique signature. 

If acoustic oscillations are present in the CMB, as is the case in all but models with significant reionization 
(see e.g. Peebles 1987, Efstathiou & Bond 1987, Coulson et al. 1994) or late formation of perturbations 
(Jaffe, Stebbins & Frieman 1995), their signature will provide information on the curvature. The damping 

tail contains the most robust information, but its location alone can constrain but not precisely fix the 
curvature. Additional restrictions such as a standard recombination history and a near BBN baryon content 
are required to make this a sensitive probe of the curvature. Furthermore, in models where the acoustic 
signature is sufficiently regular, the spacing of the peaks provides a precision test of the curvature even if 
the baryon content is anomalously low or somewhat high. It can also be combined with the damping tail to 
discriminate against truly exotic models. In §5, we discuss specifically what regularities must be observed 
before the curvature can be unambiguously measured if the correct model for structure formation is not 
assumed to be known a priori. 
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The outline of this paper is as follows. The next section contains the derivation of our principle results. 
Their impact on the features in the angular power spectrum of the CMB is described in §3. We discuss our 
main points in §§4 - 5, and present our conclusions in §6. 



2. Physical Processes 

This section develops the formalism behind our principle results and illustrates them in a series of 
concrete examples. The lessons learned here will be applied in the later sections. A summary of important 
conclusions is given in §2.9, so that the bulk of this section may be skimmed on a first reading. 

Acoustic oscillations in the CMB are inevitable if gravitational potential perturbations exist during the 
period when the Compton mean free path of a photon scattering off a free electron is much less than the 
horizon scale. In this case, the photons and electrons are tightly coupled. Since Coulomb interactions couple 
the electrons to the baryons, we refer to the system as a photon-baryon fluid. Photon pressure in the fluid 
resists gravitational compression and sets up acoustic waves in the system. 

Because the properties of the oscillator are determined by the background while those of the gravitational 
forces are described by the model for the perturbations, the acoustic signature provides a unique opportunity 
to probe both the background cosmology and the model for structure formation. Here, we explore the 
evolution of acoustic phenomena under the influence of an arbitrary source of gravitational perturbations. 
This is employed in §3 to determine the conditions under which signatures such as the peak locations, relative 
heights, and damping tail may be considered robust. These signatures will be observable in the small scale 
CMB anisotropy if the tight coupling condition is satisfied during the epoch immediately preceding the last 
scattering event. In particular, it holds for the standard thermal history where recombination and hence last 
scattering occurs at redshift ~ 10"^. 



2.1. Fluid Equations 

We start with the fundamental equations describing the dynamics of a relativistic fluid. The physical 
interpretation and description of these equations is given below. The evolution of the photons and baryons in 
a metric perturbed by density fluctuations in the fcth normal mode is given in the Newtonian representation 
as (see e.g. Mukhanov, Feldman & Brandenberger 1992, Ma & Bertschinger 1995) 

eo = -Je.-^, (1) 

Sb = -kVb - 3$, 

for the continuity equations and 

Gi = k[eo + * - in^] - f(ei - h), 

Vb = ~Vb + k^ + f (Gi - Vb)/R, 

for the momentum conservation or Euler equations of the photons and baryons respectively. Here overdots 
are derivatives with respect to conformal time rj = J dt/a, R = ip^/'ip^ is the baryon- photon momentum 
density ratio, and f = XeUeOTa is the differential Compton optical depth with as the ionization fraction, 
Ue as the electron number density, and or as the Thomson cross section. The fluctuations are defined as 
Go = AT/T = 6-^/4: the isotropic temperature perturbation, Gi the dipole moment or photon bulk velocity. 
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the photon anisotropic stress perturbation, 5b the baryon energy density perturbation and Vf, the baryon 
velocity. The gravitational sources are $, the perturbation to the spatial curvature, and 4", the Newtonian 
potential. In this gauge, these gravitational perturbations distort the metric as goo = —o?{l + and 
Qij = a?{l + 2$Q)7jj, where 7jj is the three metric on a surface of constant curvature and Q is a plane wave 
cxp(«k • x) in a flat geometry or more generally the fc-cigenfunction of the Laplacian (Wilson 1983). The 
Einstein-Poisson equations thus relate them to the matter fluctuations as (Bardeen 1980) 



pA + ^^{pi +Pi)Vi/k 



(3) 



where the sum is over particle species, the curvature K = —Hq{1 — — flj^) and the Hubble constant 
Ho = lOO/i km Mpc~^. We will assume from now on that the relevant scales are far under the curvature 

scale [K/k"^ 0). None of our main results are affected by this assumption for reasonable values of K. Note 
also that if the anisotropic stress pt^t = "^Pi^i negligible, = — 

Now let us examine the physical content of Eqs. (1) and (2). Photon number conservation relates changes 
in the temperature fluctuations to the velocity divergence by a factor of ^ (since number density is related 
to temperature by n-y T^). The $ term represents the dilation effect on the wavelength of the photons. 
Since the curvature perturbation "stretches" the spatial metric, changes in $ gives rise to a dilation effect 
of the same origin as the cosmological redshift. Notice the sign of this effect implies that the photons will 
always oppose a change in the curvature, a point which will be very important later. Similar effects govern 
the baryon continuity equation. Since the fractional energy and number density fluctuations are equal for a 
non-rclativistic particle, their rate of change is given by the velocity divergence. The additional effect due 
to the stretching of the volume from $ also implies a number density dilution of 3$. 

The expansion makes particle momenta decay as a~^. For the photons, this is accounted for by the 
temperature redshift; for the baryons, by expansion drag on the bulk velocity (the a/ a term). Gradients in the 
potential, k^, generate velocity perturbations by gravitational infall. For the photons, infall is countered by 
stress in the fluid, both isotropic (pressure) and anisotropic (quadrupole moment). The baryons however are 
effectively pressureless. The photon and baryon equations are coupled by Compton scattering (the f terms) 
which exchanges momentum between the fluids. Since the momentum density of the fluid is proportional to 
p + p, conservation relates the scattering terms by the factor R—{pi, + Pb)/{p-^ + p^) ~ ipb/'^P^- Scattering 
seeks to equalize the bulk velocities 6i = H causing adiabatic evolution of the density perturbations 
h = 390. 

If the scattering is rapid compared with the travel time across a wavelength, the momentum conservation 
Eqs. (2) may be expanded in powers of the Compton mean free path over the wavelength k/f. By eliminating 
the baryon velocity, we obtain the tight coupling approximation for the evolution of the photons (Peebles & 
Yu 1970, HSa) 

eo = -^61-$, 

The quadrupole term 11^ = 0{k/f)@i causes viscous damping and is treated in §2.8. It is a higher order 

correction because scattering tends to isotropize the photons in the baryon rest frame and suppresses the 
quadrupole. From examining Eq. (4), one can see that baryons decrease the efficacy of the pressure and add 
an expansion drag term to the momentum equation. The gravitational infall term remains unaltered since 
its baryon analogue is identical. 
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2.2. Gravitational Redshift Effects 

Before proceeding with the main task of exploring the acoustic signatures, let us review how gravitational 
effects manifest themselves in the CMB (Sachs & Wolfe 1967). As discussed in §2.1, these arc the ordinary 
redshift of a photon climbing out in and out of potential wells and the dilation effect from changes in the 
spatial metric. If the metric fluctuations are generated by the density fluctuations in the photon-baryon 
system itself, the Poisson equations (3) tell us that they are suppressed by a factor of {krj)~'^ with respect 
to the temperature fluctuations inside the horizon. Thus the seZ/-gravity of the photon-baryon fluctuations 
generally only is important outside the horizon krj <^ 1. However this is not necessarily true for metric 
fluctuations generated by an external source. Gravitational redshift effects can significantly alter the acoustic 
signature, and we must include them in the analysis even on small scales. 

Because temperature perturbations are observed only after the photons have lost energy climing out of 
potential wells ^, the "effective" temperature perturbation is given by 8o + It is this quantity that we 
measure as a temperature fluctuation on the sky. It will be important in the following sections to consider 
the effective temperature fluctuation rather than the intrinsic fluctuation @q. The blueshift from infall into 
a constant gravitational well is exactly cancelled by the redshift from climbing out. Thus 6o may have a 
large but unobservable offset which is removed in the effective temperature 6o + ^■ 

On the other hand, the intrinsic temperature evolves as Oq = — $ above the horizon to oppose changes 
in the spatial curvature. Thus the effective temperature obeys the relation 6o + = 'I' — $ ~ —2$ above 
the horizon. This yields a general description of the gravitational redshift or Sachs- Wolfe effect (Sachs & 
Wolfe 1967), 

[eo + mv,k) = [eo + mvi,k)+[^-mv,k)\l, (5) 

where rji is some initial epoch at which the fluctuations were formed. In particular, if the initial conditions are 
isocurvature ^{r]i, k) = 0, and the initial temperature fluctuation is also small, then the effective temperature 
becomes ^ — $ rj —2$. The photons are thus underdense in potential wells. For adiabatic fluctuations, 
photons are overdense in potential wells so that the effective temperature fluctuation is reduced. As we shall 
sec below, in the radiation-dominated era 8o + '5 = \^ ■ The change in the equation of state through the 
matter-radiation transition causes a small decay in the potential [see Appendix, Eq. (A4)] and brings the 
effective temperature in the matter dominated limit to 6o + ^' = (Sachs & Wolfe 1967, HSa). 



2.3. Oscillator Equation 

In HSa, a formalism was developed to calculate the response of the photon-baryon fluid to metric 
fluctuations \E' and $ which are considered external to the fluid. Combining the two equations in (4), we 
obtain 

^(1 + R)@o + yeo = -y (1 + i?)* - J-(l + i?)$. (6) 

Conceptually, this equation reads: the change in momentum of the photon baryon fluid is determined by a 
competition between the pressure restoring and the gravitational driving forces. Below the sound horizon, 

photon pressure resists gravitational compression and sets up acoustic oscillations. Though conceptually 
useful, this approach has the practical disadvantage that the photon-baryon contribution to the metric 
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fluctuations must be; already known. In general, it is unknown and we must break the metric fluctuations up 
into pieces generated by the photon-baryon fluid (76) and by the external source (s), e.g. from dark matter, 
entropy fluctuations, defects, etc. The Poisson equations (3) and the tight coupling condition VJ, = 9i tell 
us 



a 1 



Oj k 



(8) 



for adiabatic fluctuations in the photon-baryon system where 61, = SQq. As we shall see in §2.5, entropy 
fluctuations S = Sb — 30o are simply described as an external source since 5 = in the tight coupling limit. 
Thus Eq. (8) represents no loss of essential generality. 

As Eq. (8) implies, the response of the photon-baryon fluid is most easily solved in the so-called rest 
frame of the fluid (see Appendix). Here the temperature perttirbation is 



a k 



(9) 



Combining Eqs. (4) and (9), it evolves under the equations 



1 + \{1 + R) 

y 



T- 



y 



1 



y i + R 



-T 



y 



y 



6 



i--(i + i?)^ 

y2 



1 



l + R 



-T 



y ' " 

^ ft 1 



(10) 



where primes are derivatives with respect io x = krj with k fixed, and y = {fl-yHQ)~^/^ak. Notice that if 

|0i| ^ |0ol; is the case in the oscillatory regime xcg > 1, Eq. (9) implies that T = 8o[l + C>{x~^)] and 
the two representations become equivalent. In the next few sections, we will examine the implications of 
Eq. (10) through a number of examples and approximations. 



2.4. Photon Backreaction at Early Times 

The evolution equation (10) simplifies substantially in the radiation-dominated epoch. Early on, 
the baryon contribution is negligible and the fluid evolution is dominated by the response of the photon 
perturbations to the source. The resulting photon-baryon density fluctuation feeds back into the metric 
fluctuation. Even though the universe may be just becoming matter dominated at last scattering, its early 
radiation-dominated legacy plays the dominant role. As we shall see, the processes that fix the amplitude 
of the acoustic oscillation take effect mainly around horizon crossing when the universe was still radiation 
dominated for the relevant fluctuations. 

It is instructive to consider first the case in which the expansion is photon- dominated. This neglects 
the neutrino and source contribution to the background energy density, but does not fundamentally alter 
the results for the early superhorizon evolution (see HSb). In this limit, R ^ 0, y ^ x, and the evolution 
equations (10) become 

r'-r/a: + ei/3 = ^^^^ 

ei + Gi/x- [l-6/x^]T = 

or combining the two 

T = S, (12) 



T" 
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X = k77 



Fig. 1. — Pure modes. 

In the absence of sources, the growing mode of the rest frame temperature perturbation matches onto a 

cosine acoustic oscillation inside the horizon, whereas the decaying mode matches onto a sine oscillation. 
The oscillator response to an external source is constructed by Greens method from these homogeneous 
solutions and transformed into the Newtonian frame with Eq. (9). 



where the gravitational source is given by 

12 



S = 



(x2 + 6)2 



x{x^ + 18) 

+ 6)2 



x^+6 



(13) 



If 5 = 0, this is simply the Bessel equation. The homogeneous solutions are therefore (Kodama & 
Sasaki 1986) 

Ta{x) = — cos(a;/-\/3) + (a/S/x) sin(a;/-\/3), 
Tb{x) = -sin(a;/\/3) - (\/3/a;)cos(a;/\/3), 



(14) 



with Wronskian 1 / \/3. These functions are plotted in Fig. 1 and may be considered as the fundamental or 



pure modes of the photon-baryon fluid. Their limiting behavior as x is Tq = x^/9 and % 



-\/3/a 



From Eq. (9), the corresponding limits for the Newtonian fluctuations are Qq = 1/3 and 6-\/3/a;^ and 
= 2/3 and —6\/3/x^ for the two modes respectively. As a; ^ oo, they become cosine and sine waves 
respectively for both the rest frame and Newtonian temperatures. 

Although Eq. (11) may easily be solved numerically, the Green's method is more illuminating. 
Constructing the solution out of the pure modes, we find 



with 



T{x) = Aa{x)%{x) + Ab{x)%{x), 

Aa{x) = Aa{Xi)-V3j^Jb{x')S{x')dx', 

Ab{x) = Ab{xi)-VSJ^Ta{x')S{x')dx', 



(15) 



(16) 
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where Xi is the initial epoch at which the perturbations and hence the constants Aa{xi) and Ab{xi) are fixed. 
If the photon fluctuations are set to zero at Xi, then the answer can be read directly off the source behavior. 
When % ^ Ta, source contributions stimulate the 7^ mode mainly. An examination of Fig. 1 would imply 
that all superhorizon effects from the source would create a 7^ cosine mode. However, the initial fluctuations 
could be set up such that Aa{xi) and Ab{xi) exactly cancel the influence of the source. This is exactly what 
occurs in the class of isocurvature models based on balanced initial conditions. 

To see this more clearly, let us express the evolution in terms of ^-yt, the curvature perturbation 
generated by the photon-baryon fluid. In the photon dominated limit, equation Eq. (8) gives = GT/x"^ 
and Eq. (12) becomes 

xH';, + 4x^, = -xH'^ - 4x^'„ (17) 

for a; <C 1. The solution is 

Hx) = $^6(x) +$,,(.t) 

= ^-^b{Xi) + ^,{Xi) + + + {2/Xi)ei]{Xi) [{Xi/xf - 1] , 

where we have approximated ^'^(.Ti) = — $s(.Ti) in rewriting the initial conditions in terms of Oi. Aside 
from a decaying mode, the photons evolve to keep the curvature perturbation constant (c.f. Eq. (A9) and 
Veeraraghavan & Stebbins 1990). Thus, independent of the source behavior, if ^b-yixi) = —^s{xi) and 
Qi{xi) = 0, compensation forces $(a; <C 1) ~ 0. 

While this argument only shows compensation for a photon dominated system, the argument applies 
equally well for whatever the dominant dynamical component is since S cx — $ for a; ^ 0. The argument does 
not strictly apply if the identity of this component, i.e. the equation of state, changes. For example, pressure 
fluctuations and hypersurface warping can change the curvature through the matter-radiation transition even 
if a; <SC 1. We discuss these points further in the Appendix, and show that densities and hence temperature 
fluctuations arc anticorrelated with the total curvature outside the horizon for isocurvature initial conditions 
(see also Eq. 5). 



2.5. Driven Oscillations and Superhorizon Effects 

In this section, we will see how superhorizon compensation drives the oscillator and stimulates one of 
the two pure modes of §2.4. Let us start with a simple and concrete example: the baryon isocurvature case. 
Here, we begin at some initial epoch with an entropy fluctuation S{xi) = 5{nb/n-y) = Sb{xi) — 3Qo{xi). 
Tight coupling implies that the number density fluctuation of the photons and baryons, and hence S, 
remains constant [see Eq. (1)]. The entropy acts as an external source, which from Eq. (3) contributes as 
= AnGa^pbS, or rewriting this in the form of Eq. (8) 

^, = -^, = ^RS=-, (19) 

X'' X 

where A = ^{flbH^ ji^jHQ j^^^^k^^ S . Notice that the curvature perturbation implied by the source actually 
increases as x ^ 0. For this model, the effective source in Eq. (13) reduces to the simple form 

S = -l<f. = |i. (20) 

Now let us assume the isocurvature condition at the initial epoch ^{xi) = 0, or ^b-y{xi) = —^s{xi) and 
Qi{xi) = (see also Appendix). By requiring continuity in T and its first derivative, the initial partition into 
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Fig. 2. — Compensation and isocurvature fluctuations. 

(a) Baryon isocurvature model (#s oc x~^). Outside the horizon x < I, backreaction from the photons 
cancels the contribution of the source to the curvature fluctuation. Inside the horizon, pressure prevents 
significant metric contributions from the photon-baryon fluid and ^ —* ^s- (b) For source functions #s x'', 
the ratio oi Ta to Ti, amplitudes decreases due to feedback. This leaves the acoustic oscillation mainly in 
the sine mode at a; ^ 1. 

pure modes in Eq. (16) becomes Aa{xi) = —A/xi and Ab{xi) = {■\/3/54)Ax'^, yielding a large contribution 
to the Ta mode. However, let us examine the influence of the source term on the subsequent evolution, 

V3/;^ %{x')S{x')dx' = -(^/3/54)(,x2 - xj)A, 
V3j^r,ix')S{x')dx' = (x-^ -xr^)A. (x « 1) 

Thus as the evolution progresses, the initially large 7a contribution drops precipitously leaving mainly the 
Tb or sine mode at subhorizon scales a; > 1. Fundamentally, this is due to the feedback effect of Eq. (17): 
the photons oppose any change to the net curvature and evolve to maintain the isocurvature condition. 
In Fig. 2, we show the time evolution of T, the source ciu-vature and the total curvature in the baryon 
isocurvature model. By comparing $ to $s, notice that the feedback effect is only important outside the 
horizon. Another interesting quantity is the ratio of the 7^ to 7;, amplitude shown in Fig. 2b. In this case, 
the Ta mode essentially disappears after horizon crossing leaving the acoustic perturbation in a pure sine 
mode (p = —1.0, heavy line). In fact, the integral in Eq. (16) takes on a simple asymptotic form (HSb) 

lim^^oo Aa{x) = 0, 
lim^^ao Ab{x) = A/y/S 

We also show in Fig. 2b the behavior when the source is generalized to $s = —^s = Ax^. Again, the initially 
large 7^ mode is reduced as the superhorizon scale evolution progresses. 

Even if the source has some superhorizon scale feature, i.e. a maximum at .t <C 1, its effect is mainly 
cancelled out. From Eq. (18), notice that the photons always attempt to counter the source. The photon 
fluctuations would track the rise and fall of a feature leaving an effect only from the boundary conditions. 
Below the horizon however, compensation cannot occur due to the intervention of photon pressure. 



(22) 
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Fig. 3. — Driven oscillations. 

The self-gravity of the photon-baryon fluid drives a cosine oscillation for adiabatic initial conditions (thin 

lines) and a sine oscillation (thick lines) for isocurvature initial conditions. The adiabatic model has 
^s{xi) = and ^^if,{xi) ^ 0. The isocurvature model is the baryonic model of Eq. (19). The dashed 
lines show the full potential, the solid lines the effective temperature. In both cases the amplitude in ©0 + 4' 
increase during the first few oscillations, as described in §2.5. 

How then does one obtain strong contributions to the 7^ cosine mode? If superhorizon fluctuations 
are not fully compensated, Eq. (16) implies a large amplitude in 7^. This adiabatic component implies 

superhorizon curvature perturbations initially: 7^ 0. The simplest example is the case of = and 

^■yb{xi) ^ 0. In Fig. 3, this case is contrasted with the baryon isocurvature case. Let us use these examples 
to gain further intuition about the feedback mechanism. 

From Fig. 3, we see that in the isocurvature case, compensation prevents a large curvature fluctuation 
from appearing outside of the horizon regardless of the source. However photon pressure, which becomes 
important around sound horizon crossing, resists the accompanying rarefaction of the photon fluid due to 
dilation. At this point, the fluid turns around and begins falling into the potential wells of the source 
enhancing the curvature fluctuation by its self-gravity. Note the increase in the amplitude of the oscillation 
between the negative maximum and the positive maximum. As the photons resist further compression at 
the positive maximum, the self-gravity contribution to the potential fluctuation decays. This leaves the 
photon-baryon fluid in a highly compressed state and increases the amplitude of the acoustic oscillation. 
Thus the self-gravity of the photon-baryon fluid essentially drives the oscillator. It provides a kick at each 
of the first two turning points to enhance the oscillation. 

A similar analysis applies to adiabatic fiuctuations. Here the initial curvature fiuctuation is kept constant 

outside the horizon by photon backreaction. From §2.1 recall that the intrinsic temperature fluctuation and 
gravitational potential partially cancel in the effective temperature. From the a; <C 1 limit of Eq. (14), 
Go + ^■yb = ^^jb < 0. At horizon crossing, the fluid begins to compress itself due to its self-gravity and 
the effective temperature reverses sign. As pressure begins to stop the compression, the potential decays. 
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Fig. 4. — Forced oscillations. 

Inside the horizon, the regulatory effects of photon-baryon self-gravity become ineffective. A source that 

peaks at a; = fe»? 2> 1 can produce complicated forced oscillations in the fluid. The intrinsic temperature 
fluctuation T ?a ©o suffers a large unobservable offset —^s (bottom panel) which is removed by the redshift 
in the effective temperature Bo + (top panel). Notice that the fundamental period krg = 27r or i; ~ 10.9 
is still clearly apparent for this slowly-varying source. 

Again the fluid is left in a highly compressed state, and self-gravity acts as a driving term. The dilation 
due to $ (see discussion following Eq. (3)) doubles the efii'ect of infall from ^f, so that the end amplitude is 



A difFcrcncc in phase between the two modes arises since gravitational infall which initiates the chain 
of events leading to the adiabatic acoustic mode only takes effect inside the horizon, whereas the dilation 
effect that begins the isocurvature chain of events is already occurring outside the horizon. The end result is 
that the self-gravity of the photon-baryon fluid drives a cosine oscillation for adiabatic initial conditions and 
a sine mode for isocurvature initial conditions. This is not the most general case however, since self-gravity 
and thus the initial conditions only dominate the behavior around or before horizon crossing. We consider 
the behavior after horizon crossing in the next section. 



If the source is still effective inside the horizon 4>s(x > 1) > ^s{x = 1), more complicated behavior can 
result since the regulatory mechanism of photon feedback is ineffective. This is the case of so-called "active 
perturbations" (Magueijo et al. 1996). The general Greens function solution can be rewritten as 




2.6. Forced Oscillations and Subhorizon Effects 




(23) 
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for a; ^ 1. The phase of the acoustic wave is in general time dependent for subhorizon effects. Furthermore, 
we cannot ignore the gravitational redshift of the photons climbing out of the potential even inside the 
horizon (see §2.2). The effective temperature is 6o + ^' ~ T + if a; ^ 1. In Fig. 4, we plot an example 
based on a broken power law source of the form 



Here pi = 3, P2 = 1 and a = 0.01 such that the source peaks well below the horizon scale. Initial conditions 
do not play a major role here since the fluctuation is generated inside the horizon. By considering only 
the effective temperature 9o + ^E", we remove the large unobservable offset from gravitational blueshift in T 
(see §2.2). Note that it is important to include the Sachs- Wolfe contribution to obtain a reliable measure of 
the temperature anisotropy at last scattering. This is especially true when the anisotropics are "sourced" 
such that metric fluctuations remain substantial inside the horizon, as in some defect models. Neglecting 
this effect, for example by use of the oft-employed synchronous-gauge photon density, can lead to erroneous 
conclusions about the shape of the CMB anisotropy power spectrum. Should the potentials have significant 
time variability after last scattering it is also necessary to include the ISW effect, as this can lead to significant 
distortions of the peaks, notably the first peak. 

Even this smooth, well-behaved source leads to complicated structure in the acoustic oscillation. 
Furthermore, after the source turns off and Aa and Ab become constant in time, the oscillator is left with 
a phase shift described by Eq. (23) that does not necessarily correspond to one of the pure modes. It is 
important to note that even in this case, the temperature perturbation is anticorrelated with the source 
fluctuation immediately after horizon crossing. This is due to compensation from feedback around horizon 
crossing and will be important for the question of robust distinctions between models. We have found that 
this anticorrelation persists even if the source changes sign outside the horizon or near horizon crossing. 

It is of course impossible to quantify all the possibilities that arise from the arbitrary action of a source 
well inside the horizon. Howevc^r Fig. 4 illustrates two general points. Large metric and correspondingly even 
larger density fluctuations are required for subhorizon forcing to dominate over driving effects at horizon 
crossing. Compare the vertical scales of the upper and lower panels of Fig. 4. There are two effects that 
make forcing around horizon crossing easier than well inside the horizon. The first is that for a; 1 a 
large metric perturbation requires a very large density perturbation, since the Poisson equation (3) gives 
= {PsI Pt)0{6s/x^). Also, from the Euler equation, pressure perturbations are more important than 
metric perturbations at large x, once the artificial offset due to the k'i term has been removed by considering 
the effective temperature. From Eq. (6), the effective temperature evolves according to 

(Oo + *)" + i (90 + *) = - (25) 
o 

The source is thus ss — 2$", and if constant, produces a temperature shift of 6$". Exotic spectra of this 
type require extreme conditions. More plausible is a model whose forced effects dominate only around a; ~ 1 
where feedback effects still play some role. As we shall see in §3.2, this has consequences for the coherence 
of the resulting oscillations. 

The second point is that in spite of the gravitational forcing the natural period of the oscillation, here 
X = 2-\/37r is still apparent in the harmonic structure. If the source <t>" is slowly varying in time, so is the 
phase shift. In Fig. 4, the spacing between the peaks is still regular and corresponds to the natural period 
despite the slowly var}dng offset. Thus this acoustic signature is only lost if the metric fluctuations are both 
extremely large and rapidly varying inside the horizon. 
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Fig. 5. — Baryon drag and relative peak heights. 

Baryons displace the zero point of the oscillation by — ij^t. For a near BBN baryon content (center panel), 
the displacement is smaller than the oscillation itself leading to alternating peak heights in the rms. For a 
much larger baryon content, the photons oscillate around a strongly displaced zero point. 

If such an exotic spectrum is measured in the CMB, we would have concrete evidence that the mechanism 
responsible for structure formation requires new physics. Until the spectrum is measured however, we can 
only consider the broader implications. The possibility of driven sub-horizon effects raises two questions 
that we will address in §§4 - 5. Are there any unique signatures of the pure cases that distinguish them 
from these more exotic scenarios? How does the existence of such exotic models affect our ability to measure 
quantities such as the curvature of the universe? To help answer these questions, we need to understand two 
additional acoustic effects that can provide clues to unraveling the spectrum. 

2.7. Baryon Drag 

Up until this point, we have neglected the effect of the baryons in the photon-baryon fluid. This 
is appropriate for the early evolution and reveals the qualitative structure of the acoustic signal even at 
recombination since i?(?7*) = 3/9t/4p-y|^, w 31.5riiJi^ ^ 30%. Nevertheless, there is an important acoustic 
signature associated with the baryons if Qbh^ is near or greater than this BBN value (HSa, and Hu & 
Sugiyama 1996, hereafter HSc). To include its effects properly, we must solve the full Eq. (10). However it 
is instructive to examine the qualitative origin of the effect first. Eq. (6) tells us that the baryons contribute 
to the incrtial and gravitational mass of the fluid. Thus with a higher baryon content, gravity can compress 
the fluid more strongly inside the potential well 4* < 0. In the limit that a; ^ 1 and we observe for only a 
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time short compared with changes in R and the source, the solution to the oscillator equation (6) is 

Oo + * = Ci cos{krs) + C2 sin(fcrs) - i?*. (26) 

The last term represents the baryon drag on the photons and displaces the zero point of the oscillation. 

Since — > inside a potential well, it enhances the compression and suppresses the rarefaction stage of 
the acoustic oscillation in the potential well. The crucial point is that the presence of baryons allows us to 
distinguish between the two stages from the observable rms deviation |9o + Fig. 5 shows an example 
based on the source in Eq. (24) with pi = 2, p2 = 1 and a = 0.05. Here we take the A;-spectrum at 
recombination for a baryon content of Q^bh^ ~ 0, 0.0125, 0.25 where i?* = R{ri^,) = 0, 0.38, 7.6. The neutrino 
and CDM densities have been set to zero for simplicity. Notice that if baryon drag is significant but not 
dominant (middle panel) then it will modulate the peaks into a pattern of alternating peak heights. If it 
dominates (bottom panel), the rms fluctuations no longer possess zero crossings but oscillate around some 
large d.c. offset. Either pattern is quite distinctive and as we shall see can separate the adiabatic from 
isocurvature scenarios in all cases save those with <C |6o + ^I/l at last scattering. 



2.8. Photon Diffusion 

Finally, photon diffusion leaves a robust signature by providing a cut off scale to the acoustic oscillations 
that is independent of the source fluctuation. As the diffusion length passes the wavelength, acoustic 
oscillations are exponentially damped (see Silk 1968, Weinberg 1972). Physically this occurs since diffusion 
erases temperature differences across a wavelength and causes viscosity (or anisotropic stress) in the fluid. 
Anisotropic stress, or the quadrupole, is generated from the free streaming of a dipole fluctuation Oi. As 
photons from crests and troughs of the original velocity perturbation meet, their Doppler shifts create 
a quadrupole temperature pattern. This transfer of power from the dipole to the quadrupole is but a 
manifestation of a general tendency. Streaming transfers power to higher angular moments since the original 
temperature fluctuation subtends a smaller and smaller angle as seen by a distant observer. However in 
the tight coupling limit, streaming is coUisionally suppressed by the factor r/fc (the optical depth through a 
wavelength) such that in the tight coupling limit, the anisotropic stress is approximately (see HSc, Eq. A8) 

Il^ = l{k/f)f^'eu (27) 

where /2 = | (Kaiser 1983). Other approximation commonly employed are /2 = ^ for unpolarized 
radiation (Chibisov 1972) and /2 = 1 for further neglecting the angular dependence of Compton scattering 
(Weinberg 1972, HSa). Prom Eq. (26), it is natural to try a solution of the form expi Jcodrj for both 
[Go + (1 + R)'^] and 61. Heat conduction, proportional to Vb — 81, is described by iterating the baryon 
Euler equation (2) to second order, 

Vb-ei = -f-'^R[iojei - fc*] - f-^iJ^w^Oi. (28) 

Here we have ignored changes on the order of an expansion time compared with those at the oscillation 
period. Combining Eqs. (27) and (28) in the photon Euler equation (2), we obtain the dispersion relation 



1 2 -1 
±kcs + i-k f 
6 



R 4 1 1 



(l + i?)2 + 5'^2 



(29) 



where recall that the sound speed Cg = 1/ y'3(l + R). It follows that acoustic oscillations are damped as 
(HSc) 

60 + * = [60 + ^']e-[*^/'=°('')l' + i?^f(e-l'=/'=°('')l' - 1), (30) 
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where 9o is the acoustic signal described by Eq. (4) and the damping wavenumber is 

To order of magnitude, the diffusion length is the geometric mean of the Compton mean free path and 
the horizon rj, as expected of a random walk. Aside from the residual baryon drag effect, photon diffusion 
leaves a distinctive damping signal in the CMB that is only dependent on the cosmological background. 
Because the damping is exponential, no oscillations survive at ko regardless of the source. 

One complication arises however. As recombination progresses, the ionization fraction and hence 

the differential optical depth f decreases. The corresponding increase in the damping length can still be 
approximated in the tight coupling limit if f/kDiv*) ~ kar]* ^ 1. In this case, the diffusion length passes 
the wavelength while tight coupling holds and the damping can be calculated semi-analytically from the 
known ionization history (see HSc). If this condition is not satisfied, as is the case for extremely low ilfjft,^ 
or reionized models with a long Compton mean free path at recombination, the damping must be calibrated 
numerically (Hu & White, in preparation). 



2.9. Summary 

In summary then, we have shown that the evolution of fluctuations acts to resist changes in the spatial 
curvature, whenever the pressure can be neglected (outside the horizon). Schematically, this is because 
the radiation or matter in a growing curvature perturbation redshifts faster than average (c.f. cosmological 
redshift, as in §2.1) reducing the local energy density and hence stabilizing the curvature perturbation. It is 
also a consequence of causality as shown in the Appendix. 

The behavior of the acoustic oscillations due to sources outside the horizon thus depends only on 
the initial conditions: whether or not there are initial uncompensated curvature perturbations. The 
compensation by the radiation also implies that the self-gravity of the photons produces an important 
feedback mechanism. The result is that sources arising before or near horizon crossing stimulate either the 
cosine or the sine mode of the oscillation, not an arbitrary admixture. 

The fact that the compensation mechanism leads to rarefactions in the photon density at horizon 
crossing for isocurvature fluctuations but compressions for adiabatic fluctuations inside the potential well is 
an extremely robust conclusion. Of crucial importance is it leads to an observable effect in models with f2(,/i^ 
of a few percent or greater due to the drag induced by the baryons. Photon diffusion also leaves a distinct 
damping tail that is entirely independent of the model for the fluctuations but sensitive to the thermal 
history. 



3. Acoustic Signatures 

Having studied the acoustic oscillations in the last section, let us now consider what signatures the 
effects summarized in §2.9 leave in the angular power spectrum of CMB anisotropies. Acoustic oscillations 

are frozen in at last scattering r]» when the Compton optical depth to the present drops below unity. The 
evolutionary properties of the acoustic phenomena before last scattering influence the peak locations, heights, 
and damping tail. These will be used in the next sections to devise tests of the model for structure formation 
and the background cosmology. 
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3.1. Peak Locations for Driven Oscillations 

At last scattering, each fc-mode is caught at a different phase in its oscillation. Of course, one can adjust 
the magnitudes in k by choosing different relative weights for the fc-modes at the initial conditions. We 
assume that any such weighting is not in itself oscillatory in k. Under this single assumption, the acoustic 
pattern of peaks is robust to the details of the model even though their absolute heights are not. These 
oscillations on the last scattering surface are viewed today as anisotropics on the sky. Thus the fc-space 
power spectrum is projected onto an angular power spectrum. The angular correlation function is broken 
up into £ Legendre moments which represent the power on the scale 9 ^ The power spectrum is usually 
denoted as £{£+ l)Ce where is the ensemble average of the squared multipole coefficients. 

As we shall see in §5 [(see Eq. (39)], since the projection of k onto i is highly dependent on the curvature, 
knowledge of the physical scale of acoustic features allows a sensitive probe of the curvature. The physical 
scale of the peaks is related to the sound horizon = J Cgdr] ^/\/3 as — > or more generally 

r.--^(f2oi?o) ^^^In ^^^^= , (32) 

assuming the universe was radiation dominated before a^q. The exact relation of the peak scale to the sound 
horizon depends on the nature of the fluctuations and supplies a test of the model (see §5.6). Let us now 
examine this relationship more carefully. 

Once the source becomes ineffective |0o| 3> \^s\, the partition into pure mode amplitudes Aa and At 
of §2.4 becomes time independent and the acoustic oscillations settle into a pure sinusoidal form, sin(/)(?7, k), 
with phase 

(t) = kr,+ tan-^{Aa/Ab). (33) 

For sources which peak before horizon crossing, backreaction effects create a two state system (see §2.5). If 
the fluctuations are not perfectly compensated at the initial conditions, Aa '> Ab- Compensation tends to 
create a situation where At > Aa- 

Robust predictions arise for the location of the peaks in the k or i space power spectra. Since the phase 
is dependent only on the nature of the initial fluctuation, not on the detailed behavior or the source in time 
or k, all such models give definite predictions for the peak locations. 

(f){k,„) = {m-l/2)TT, (34) 

for integer m. Furthermore, the harmonic series of peak locations is independent of the background quantities 
that fix the sound horizon and the angle it subtends on the sky. The generic adiabatic prediction is that peak 
locations follow a 1 : 2 : 3 • • • series in k or i. Isocurvature models tend toward a 1 : 3 : 5 • • • pattern (HSb). 

Even if Aa and At are comparable, as would be the case if a small (finely tuned) uncompensated fluctuation 
remained in the initial condition, the pattern of peaks uniquely determines the phase shift. This might be the 
case if the initial conditions contain a balanced amount of coherent adiabatic and isocurvature fluctuations. 
In most physical examples however (see e.g. Kawasaki, Sugiyama & Yanagida 1995), the processes which 
generate the two types of fluctuations are statistically independent and the two contributions are incoherent, 
i.e. generated with different phases in k. In this case it is even more unlikely that the two contributions 
would be balanced to give the same acoustic amplitudes. For incoherent superpositions, the phase is always 
determined by the dominant component with a change in amplitude, not phase, caused by the weaker 
component. Furthermore, the phase difference <j){km) — (t>{km-i) = implies that regardless of the phase 
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shift or amplitude variation, the spacing of the peaks in k is given by 

kA = km- km-l = T^/rs (35) 

Thus the acoustic peaks possess both model information in the ratios of their locations and model independent 
information for the measurement of background parameters in the spacing between the peaks. 

3.2. Peak Locations for Forced Oscillations 

Although dcfinitcncss in phase is a typical feature of acoustic oscillations, it is not necessarily obeyed by 
models where the source is still active inside the horizon. In this case, the partition into pure modes through 
Aa and At is time dependent. There are two limiting cases worth considering. The source could be scale 
free in its temporal behavior such that k) = G{x)F(k) where recall x = kr]. This occurs for example if 

the source behavior is correlated with horizon crossing due perhaps to the onset of a causal mechanism. As 
discussed in §2.6, this is in fact the most likely scenario since metric perturbations that can overcome photon 
pressure become increasingly hard to generate inside the horizon. In this case, the behavior as a function 
oi X = krj represents both the time evolution of a fixed fc-modc and the transfer function in k for a fixed 
time, e.g. at last scattering 77*. If G{x) peaks around x = 1 as is typical, we would expect an irregular first 
peak followed by an increasingly regular but phase shifted harmonic series. If it peaks much after x = 1, 
complicated acoustic behavior containing detailed information about the source evolution would result (see 
Fig. 4). Even in this case, we can see from Eq. (25) that if the source is slowly varying in time such that 
G"" <C c^G", the natural period is still imprinted in the oscillations. In cases where the driving force is 
large, we expect only compressional phases to be clearly visible as peaks in the rms. Rarefactions become 
rms troughs. Using only the compressional extrema, the peak spacing is then given by 2kA = 27r/rs. 

In the extreme case that the metric perturbations are generated well after x = 1, one might also expect 

stochastic behavior in the source (Albrccht ct al. 1995). Here the timing mechanism due to the act of horizon 
crossing no longer serves to correlate the modes. In this case, each fc-mode receives a different set of impulses 
from the source. The phase shift given by tan~^[^o(a;, A:)/^f,(a;, A:)] varies with k. In the extreme limit of 
rapid variation, the phase information is lost when summing over fc-modes to form the observable anisotropy. 
This washes out the oscillatory behavior. It is believed that a concrete example of this mechanism is given 
by the cosmic string scenario (Albrecht et al. 1995) although whether it is the dominant mechanism or not 
is currently disputed. We shall discuss what can be learned in the case where the oscillations are washed 
out in §5. Let us assume for the present that oscillations will be observed. 

3.3. Peak Heights 

A wealth of information is stored in the peak heights. Their signature is more model dependent than 
the peak locations and provides an excellent means of examining the fine details of the model. On the other 
hand, robust features that distinguish between general classes of models are more difficult to isolate. In 
this section, we will examine generic features in the peak heights that may at least be used as clues for this 
purpose. 

Let us begin by examining the familiar scale-invariant adiabatic case. In the tight coupling limit, there 
are two effects that determine the heights of the peaks: the driving force of feedback and baryon drag. 
They are to a certain extent mutually exclusive. As we have seen in §2.5, feedback boosts the oscillation 
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amplitude above the initial conditions [60 + \E'](77i, fc) = ^^{r]i,k) by a factor of 3. Recall that the fluid 
first is compressed by its own self-gravity. Photon pressure resists the compression, causing the photon- 
baryon contribution to the potential to decay. The fluid is thus released exactly in this highly compressed 
state into the acoustic phase. We show in the Appendix that the potential in the matter dominated limit 
^'(?7, k) = ^^'(774, k), thus the boost represents a factor of 5 enhancement over the 90 + ^* = Sachs-Wolfe 
tail. The inclusion of the neutrino background slightly lowers the amplitude (see HSc). The driving effect 
does not occur if the potentials are dominated by an external source such as cold dark matter. Hence the 
prominence of the acoustic oscillations increases if the universe is made more photon-baryon or radiation 
dominated at last scattering. 

Baryon drag enhances the compressional, here odd, peaks by a term of 0(7?^). This alters the peak 
heights to give the distinctive alternating or offset oscillation pattern but becomes sub-dominant if the 
potential is intrinsically small. Such is the case for driven oscillations since the potential ^ w and 
decays inside the sound horizon. 

For isocurvature models, gravitational redshifts cause the effective temperature above the horizon to 
be 00 + '5 » 2'!' < (sec §2.2). ft is also worthwhile to note that although the fluctuation outside the 
horizon at last scattering is small in these models, this does not imply that the observable anisotropy from 
those scales is correspondingly small. Gravitational redshifts from the time dependent potential continue to 
generate fluctuations between last scattering and the present in the same manner. This is generally called the 
integrated Sachs- Wolfe (ISW) effect. Thus an isocurvature model which generates a scale invariant spectrum 
of curvature fluctuations near horizon crossing yields a flat large scale anisotropy spectrum just as in the 
familiar adiabatic case (Bennett, Stebbins & Bouchet 1992, Coulson et al. 1994, Pen, Spergel & Turok 1994). 

Superhorizon isocurvature evolution makes the photon-baryon fluid more and more rarefied inside 
a potential well until photon pressure can successfully resist rarefication. Since the first feature is the 

turning point from the superhorizon behavior, it is not prominent in comparison to the 2^* Sachs- Wolfe 
effect. However, the fiuid then begins infall into the source wells. The driving effect of the photon-baryon 
contribution to the potential wells now proceeds as in the adiabatic case to enhance the oscillation making 
the second feature much more prominent than the first (sec Fig. 3). Whereas the first peak has a height of 
order 25' at kvs — tt/2 and is small since compensation eliminates metric fluctuations above the horizon, the 
second peak has a height of order 2\1/ at krs = 37r/2, which is significantly larger since the photon-baryon 
contribution adds to rather than cancels the source. As in the adiabatic case, baryon drag contributes an 
0{R^) term that boosts the compressional phases. However in the isocurvature case, these are the even 
peaks since compensation demands that if the first feature occurs near the horizon it represents rarefaction 
inside the potential well. Thus, the second isocurvature peak is lifted even higher with respect to the first 
by baryon drag. 

How robust are these general tendencies? The prominence of the acoustic oscillations compared with the 

large scale tail can be masked or altered by the presence of tilt or features in the initial /c-spectrum as well as 
by other effects from tensor and vector metric perturbations (e.g. Crittenden & Turok 1995). These effects 
are however unlikely to obscure the distinctive alternating peak heights due to baryon drag. Unfortunately, 
the baryon drag effect may be difficult to observe if \R^\{r]*) <^ |6o + ^Hv*) ~ l^liv ~ l/cgk). Even with 
high precision measurements, one must first remove the diffusion damping envelope in this case. 
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Fig. 6. — Diffusion damping. 

Although adiabatic and isocurvature models predict acoustic oscillations in different positions, they both 
suffer diffusion damping in the same way. The damping length is fixed by background assumptions, here 
f!o = 1, /i = 0.5, r^i, = 0.05 and standard recombination. These calculations were performed using a 
full numerical integration of the Boltzmann equation with the code of Sugiyama (1995) as were results in 
Figs. 7,8,10,11. 

3.4. Damping Tail 

Unless baryon drag dominates over the acoustic oscillation, as in the case that |i?^'|(77*) 3> |6o + (??*), 
the damping leaves a clear signal in the CMB. To demonstrate this robustness to changes in the model for 
the gravitational source, we compare in Fig. 6 the anisotropy power spectra of a standard inflationary model 
with scale invariant curvature fluctuations at horizon crossing and a similar axionic isocurvature model. The 
background parameters are set to be equal at = l.Q, h = 0.5 and Q.b = 0.05. Notice that the damping 
behavior is independent of the nature of the fluctuations. 

The angular location of the damping tail is highly sensitive to the curvature which dominates projection 
effects (see §5.3) but also to the thermal history which sets the maximum diffusion scale, the baryon content 
Vthh? which sets the mean free path of the photons, the matter content Vloh"^ which sets the horizon scale. A 
measurement of the damping scale alone would fix the combination of these quantities defined by Eq. (31). 
More specifically, if last scattering occurs due to standard recombination, the wavelength A;^ at which the 
acoustic amplitude falls to of its original amplitude can be fit by semianalytic techniques (HSc, Eq. E6) 
to ~ 10%. Unfortunately, a measurement of ku cannot break the degeneracy in these quantities. If Q,},h^ 
is constrained by BBN and ^qK^ by dynamical mass and Hubble constant measurements, it would provide 
interesting constraints on the curvature of the universe as we shall show in §5.3. Furthermore with a 
measurement of the acoustic peaks themselves, the degeneracy can potentially be removed in many models 
for the source fluctuations (see §5.6). 
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4. Uniqueness of the Inflationary Spectrum 

The inflationary paradigm is the front running candidate for a mechanism of fluctuation generation in the 
early universe: the perturbations in density which are the precursors of galaxies and CMB anisotropics today. 
Are there unique signatures of inflation that can validate the paradigm? For a long time it was thought that 
a nearly scale- invariant spectrum of anisotropics would provide evidence for inflation (Kolb & Turner 1990). 
But then it was realized that models of structure formation based on topological defects also naturally formed 
scale-invariant spectra (see §3.3 and Bennett, Stebbins & Bouchet 1992, Coulson et al. 1994, Pen, Spergel & 
Turok 1994). From the other direction, finding that the universe had non- vanishing spatial curvature once 
seemed like a way to "disprove" inflation, however recently it has been shown that inflation could survive 
such a revelation (Bucher, Goldhaber & Turok 1994, Yamamoto, Sasaki, & Tanaka 1995). 

Currently the most popular means of "proving" inflation is to test the consistency relation between 
the tensor (gravity wave) and scalar (density) modes. By measuring the detailed shape of the anisotropy 
spectrum, one may infer the relative amplitudes of scalar and tensor perturbations. The ratio of tensor 
to scalar contributions is proportional in slow-roll inflation to the slope of the tensor spectrum (Liddle & 
Lyth 1993, Davis et al. 1992, Turner & White 1996). However this method of proof requires a precise 
determination of the tensor to scalar ratio, which in turn requires a sizable fraction of the anisotropy be 
contributed from tensor modes (White, Krauss & Silk 1995, Knox & Turner 1994). What one really desires 
is a test based on the most basic ideas of the inflationary scenario that will be both observationally feasible 
and able to survive the ingenuity of model builders. 

The key feature of inflation for our purposes is that it provides a mechanism of connecting parts of 

the universe at early times which are currently space-like separated. In fact, it can be shown that inflation 
is the unique causal mechanism for correlating curvature perturbations on scales larger than the horizon 
(Liddle 1995, Hu, Turner & Weinberg 1994). The possibility of a white noise spectrum of superhorizon 
curvature perturbations is discussed in the Appendix, but is already observationally ruled out as a means 
of structure formation. The question therefore arises, are there unique consequences of such super-horizon 
curvature perturbations? If so, their observation would provide a "proof" of the inflationary paradigm. 



4.1. Robustness of the Harmonic Series 

We have seen in §3.1 that the harmonic series of acoustic peaks in a model with super-horizon curvature 
fluctuations is given by the cosine scries 1 : 2 : 3 • • • for the locations of the peaks. Compensated super- 
horizon fluctuations follow the sine series 1 : 3 : 5 • • •. Thus by measuring the ratio of the first three peak 
positions in i space, these two possibilities may be distinguished. There are two concerns that need to be 
addressed for this potential test of inflation. How robust is the harmonic prediction in the general class of 
inflationary models? Can an isocurvature scenario where fluctuations are generated inside the horizon mimic 
an inflationary series? 

Let us consider the first question. Residual driving effects can distort the shapes and locations of the 
first few peaks. Even for the pure adiabatic mode of §2.4, the peaks have not completely settled into their 
asymptotic forms until krg '> 1 and the peak positions follow the series in 0o + 

0.88 : 1.89 : 2.93 : 3.94 : 4.96 • • • . (36) 

We will hereafter refer to these ratios as the "cosine series" despite the fact that it has not yet converged 
upon an actual cosine oscillation by the first few peaks. In particular, the first peak is somewhat low in 
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Fig. 7. — Relative peak locations and the harmonic series. 

In the ideal case, inflationary acoustic oscillations follow the cosine series and driven isocurvature models 
a sine series (light dotted lines). The ISW effect, baryon drag and diffusion damping serve to distort the 
peak locations. The isocurvature cases considered are: baryon isocurvature (HSb); textures (Crittenden & 
Turok 1995); axionic isocurvature (Kawasaki, Sugiyama & Yanagida 1995) and hot dark matter isocurvature 
(de Laix & Scherrer 1995). Numerical calculations (points) are normalized to the ideal predictions of 
Eqs. (37) and (38) at the third peals and are specifically for + f^A = 1 though this constraint is irrelevant 
for the peak rations. They demonstrate that the two cases remain quite distinct especially in the ratio 
between the first and third peak. 



£ compared with the expectation from the higher peaks. Let us compare this prediction against actual 
inflationary models from a numerical calculation (see Fig. 7, left panel, solid lines). All models do indeed 
exhibit the harmonic structure predicted by Eq. (36). There is a slight shift of the first few peaks downward 
in £ as floh^ is lowered. As noted by HSb, if the radiation is still dynamically important after last scattering, 
the resultant ISW effect will shift power toward larger angles and distort the first peak. The even peaks are 
also shifted upwards as flgh^ increases to make potentials deeper and baryon drag more important. Still, 
for reasonable matter content and a near BBN baryon content a ratio of 3rd to 1st peak of £3/^1 ~ 3.3 — 3.7 
is a robust prediction of inflation. Furthermore the ratio of the first peak location to the peak spacing is 
between £i/£a ~ 0.7 — 0.9. The peak spacing is best measured by the prominent peaks, e.g. for inflation 
£a = {£3 — £i)/2. The only caveat is that models with extremely high ^Ibh"^ > 0.04 and flo ^ fi;, may be 
dominated by the baryon drag effect. In this case, only the compression phases are visible as peaks in the 
rms. Thus every other oscillation is unobservable and the series becomes 1 : 3 : 5 • • • like the isocurvature 
spectrum. For a model with reasonable matter content floh^ < 0.25, no value of ilbh"^ entirely eliminates 
the second peak. For HqH^ = 0.64, the second peak disappears for 0.04 < Ofc/i^ < 0.5 and the fourth for 
0.128 < rifeft,^ < 0.32. Even in this extreme case, only a gross violation of BBN would entirely mask the 
inflationary pattern. Thus the harmonic series is a robust prediction of inflation. 
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Fig. 8. — Obscured isocurvature peak. 

The first isocurvature "peak" appears as a shoulder and may be obscured especially in high baryon models 

where the second peak is significantly higher (left panel, arbitrary normalization). If the second peak starts 
the harmonic series, the ratios (right panel, points) can be quite close to the cosine prediction (dotted lines). 
The points are normalized to the cosine prediction at the third peak. These models can be distinguished by 
the peak to spacing ratio and the morphology of the first compressional peak (2nd feature) . 



4.2. Uniqueness of the Harmonic Series 

Is the cosine harmonic series a unique prediction of inflation? Just as in the inflationary case, for 
isocurvature models residual driving effects create a downward shift from the 1 ; 3 ; 5 • • • pattern. Even a 
pure driven isocurvature model with a source oc x~^, which converges to the sine series for fcr^ 1 (see 
§2.5), gives 

0.85 : 2.76 : 4.83 : 6.89 : 8.91 ••• . (37) 

We will hereafter refer to these ratios as the "sine series" . How closely do real models follow this prediction? 
For comparison, we show in the right panel of Fig. 7 the ratios of the peak locations for common isocurvature 
models found in the literature: the baryon isocurvature (HSb), textures (Crittenden & Turok 1995), axionic 
isocurvature (Kawasaki, Sugiyama & Yanagida 1995), HDM isocurvature (dc Laix & Scherrer 1995). Notice 
that these models roughly correspond to the sine series prediction. Thus in all inflationary models with 
reasonable baryon content, the ratio of peak locations should be distinguishable from the current models 
based on the isocurvature restriction. This suggests that to mimic inflation an isocurvature model would 
have to either possess extreme conditions or fine tuning. 

Let us attempt to construct models which mimic inflation. Since the first inflationary peak appears 
immediately after sound horizon crossing fc = Tr/r^, it is difficult to avoid or counter the driving effects which 
produce the sine series. Can a series based on peaks arising after sound horizon crossing somehow mimic the 
cosine series? Consider the possibility that the first isocurvature peak in the sine series Eq. (37) is missing 
or cannot be observed. This is indeed likely in some models since the first peak can be quite low in intrinsic 
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amplitude (see Fig. 8, §3.3 & §4.3). Furthermore, external effects such as tensor and vector modes which 
naturally possess a feature on the horizon scale at last scattering might mask or distort the peak. In this 
case, the observable peaks would follow the series: 2.8 : 4.8 : 6.9 • • • which is close enough to the cosine series 
to cause some concern, especially if only the first two peaks are measured. More generally, one might have 
a source that turns on only after sound horizon crossing so that the harmonic scries is again shifted toward 
smaller scales. These possibilities however create spectra which are not as close to the inflationary prediction 
as they initially might appear. 

The crucial point is that although the starting point of the harmonic series can change with the model, 
its spacing cannot. The separation between the peaks is fixed by the sound horizon at last scattering. 
Thus the harmonics cannot simply be scaled to match the inflationary prediction. The distinction is clearer 
when we consider the ratio of the peak location to peak separation. The idealized inflationary prediction of 
Eq. (36) requires a ratio of the first peak to the separation between peaks of 1\IIa ~ 0.88. If the second 
peak of the isocurvature prediction Eqn. (37) is to be taken for a first peak, it gives a corresponding ratio of 
('i/^A ~ 1.33. In other words, we expect a factor of 1.5 difference between the two cases. The uncertainty 
displayed in the real models of Fig. 7 do not destroy the test. When comparing models with the same 
values for Q.Qh? the idealized expectation is validated, e.g. between infiation and axionic isocurvature with 
Q.oh? = 0.25 and l^bft.^ = 0.0125 there is a 1.4-1.5 difference in ii/iA depending on which peaks are used to 
measure the spacing. Even if the background parameters are unknown, to get the ratio as low as it is in the 
highest inflationary case [IxjiA ~ 0.9 for ^qY? < 1), the first peak must be created around sound horizon 
crossing. 

The robustness of this test arises because the natural frequency of the oscillator is related to the sound 
horizon. Thus the oscillations, if regular enough to be mistaken for inflation, must appear regularly spaced 
with respect to that scale. It is of course in principle possible that the external source is periodic and 
drives the oscillator at a different frequency. It is amusing to note that even this rather unlikely scenario 
can be distinguished. The damping tail provides a scale that is entirely independent on the nature of the 
fluctuations. If the background parameters are known, the ratio of the damping scale to peak position 
provides an additional consistency test for the inflationary scenario. In models where the fluctuation is 
generated inside the sound horizon, the ratio ixjli) increases. In particular, for the above case of a missing 
first isocurvature peak, the ratio increases by a factor of 1.5 exactly as with the case of the peak-to-spacing 
ratio. 

Therefore, the only way to mimic the inflationary series is to tune the behavior of the source at sound 
horizon crossing so that it immediately generates a peak in the proper position but leaves no residual effects 
that would distort the series of higher peaks. Although contrived, this is possible. In Fig. 9, we plot a source 
given by Eq. (24) with pi = 2, p2 = 1, and a = 0.05 with f2(, — > 0. In this case, we have constructed a source 
which dies off after x w oT^ such as to leave the oscillations in a nearly pure cosine mode. The first five 
peaks follow a series 1 : 2.05 : 3.16 : 4.23 : 5.29, very close to the canonical infiationary prediction. Notice 
however that the prediction is 180 degrees out of phase with the inflationary prediction. This is an important 
fact that we will make use of in the next two sections. As we shall see, the crucial point is the effect of 
compensation near sound horizon crossing. One should also bear in mind that we have tuned the source 
to produce as pathological a case as possible. The cosine harmonic of peaks is essentially but not entirely 
unique to the inflationary paradigm. Its confirmation would strongly support the inflationary scenario. 
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Fig. 9. — Pathological isocurvature model. 

Here forced isocurvature fluctuations are tuned to match the locations of the inflationary prediction (upper 
panel dotted lines) with vanishing baryon content R — > 0. Notice that even in this case the isocurvature 
oscillations are out of phase with the inflationary prediction by 180 degrees. With the inclusion of baryon 
drag, this leaves an observable signal in the rms. 



4.3. First Peak 

Can we employ additional information to eliminate the possibility that an isocurvature model might 

mimic the inflationary model through the peak positions? Compensation from photon feedback discussed 
in §2.4 and §3.3 provides an essential distinction. Near or above the horizon, the photons act to resist 
any change in curvature produced by the source if the universe is radiation dominated (see Appendix for 
relaxations of this assumption). Thus, the first peak in an isocurvature model, if it is sufficiently close to the 
horizon to be confused with the inflationary prediction, must be anti-correlated with the source. In other 
words, the first peak in the rms temperature represents the rarefaction stage inside the potential well of the 
source rather than a compression phase as in the inflationary prediction. 

Since the dilation effect Oq = — which causes 60 + ^' ~ —2$, creates both the first peak and the 
Sachs- Wolfe effect, the first feature in an isocurvature spectrum is not truly a peak but a smooth turnover 
from the large scale behavior as discussed in §2.5 and §3.3. This point represents the epoch at which the 
fluctuation ceases to follow the growth of the curvature and turns around to start the acoustic oscillation. 
It creates the shoulder appearance of the flrst isocurvature feature (see e.g. Fig. 6 and Fig. 8). On the 
other hand, the inflationary fluctuation in the effective temperattne is proportional to cos(fcrs) and passes 
through a zero at kr^ = 7r/2 before the first peak (sec Fig. 3). The resultant spectrum thus exhibits a rather 
sharp break between the gravitationally dominated Sachs- Wolfe tail and the first peak (see Fig. 6). The 
gravitational nature of the first peak in an isocurvature model makes it substantially less prominent when 
compared with the low-£ tail. The first isocurvature peak also tends to be low with respect to the higher 
peaks because it appears too close to horizon crossing to experience the full forcing effect by the source. 
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Fig. 10. — Baryon drag in the inflationary model. 

Baryon drag enhances the compressional, here odd, acoustic peaks, (a) Although diffusion damping at small 
scales coupled with the intrinsically small baryon drag effect in low f2(,fe^ models hides the effect, the third 
peak is clearly anomalously high in all but the most extreme case Qf^h'^ = 0.0025 which is in clear violation 
of BBN constraints, (b) Lowering Qq'*^ also reduces the effect by reducing the potential fluctuation ^. 

Unfortunately these morphological distinctions are hard to quantify. Indeed, the sharp rise to the first 
peak in the inflationary model can be masked by the presence of the integrated Sachs- Wolfe effect in a low 
f2o/i^ universe (HSa). In an isocurvature model, the prominence of the first feature can be enhanced by 
a smooth bend in the power spectrum or vector and tensor modes. Nevertheless, the physical distinction 
between the first adiabatic and isocurvature peaks does suggest a robust way of distinguishing the models 
as we shall now discuss. 



4.4. Relative Peak Heights 

Whereas the inflationary spectrum obeys a compression-rarefaction-compression pattern an isocurvature 
model displays a rarefaction-compression-rarefaction pattern. In the absence of baryons, there is no 
observable distinction between compression and rarefaction since only the rms can be measured. However, as 
we have seen in §2.7, baryons enhance the compression at the expense of rarefaction leading to an alternating 
series of peaks in the rms. In fact, the example in Fig. 5 is the same pathological model of §4.2 which has 
peaks at the inflationary locations by construction. Notice that for reasonable baryon content, the even 
peaks are enhanced by the baryons, whereas the odd peaks are enhanced under the inflationary paradigm. 
Since this non-monotonic peak pattern cannot be reproduced without introducing the appropriate features 
in the initial fc-spectrum of fluctuations, the pattern of anomalously high odd cosine peaks is a unique feature 
of standard inflationary models. 

Of course, an isocurvature model can also begin with a compression if the first peak is so low in amplitude 

as to be unobservable (see Fig. 8). However, such a model cannot simultaneously create the peak positions 
of the inflationary prediction. In fact, regardless of the peak positions, isocurvature models are still unlikely 
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to mimic a strongly alternating heights pattern. Since the first compressional isocurvature peak occurs well 
after horizon crossing, neither the potentials of the source nor the photon-baryon backreaction are large 
enough to cause significant baryon drag unless the baryon fraction is exceedingly high. In this case, the 
behavior of the potential implies an anomalously high first compressional peak compared with the second 
due to baryon drag from the residual self-gravity (see Fig. 8). 

Although this pattern of peak locations and relative heights is unique, it is not an entirely robust 
prediction of inflation. The intrinsic fractional effect at the mth peak is on the order 

^(^*)|r^'^ « R{rj^)T{mkA) iv* > »7eq) (38) 

where T{k) is the matter transfer function (Bardeen et al. 1987). However, if the intrinsic effect is small, 
complications such as possible smooth tilts in the initial spectrum and diffusion damping will make the signal 
difficult to observe (see Fig. 10). For fio^^ ~ 0.25 this occurs if fifc/i^ < 0.007 and for BBN baryon content 
for flcih^ < 0.1. Still, with high precision experiments, a sufficiently smooth initial spectrum, and known 
thermal history, it is possible that even these extreme cases may be testable (Hu & White, in preparation). 
On the other hand if R{r]^)T{mkA) '> 1, baryon drag dominates and the even m peaks (rarefactions) no 
longer appear as peaks in the rms (see Fig. 5). In this case, the oscillations are small compared with the 
offset and the distinction between compression and rarefaction stage is that they are maxima and minima of 
the rms respectively. Since T{k) < T(0), a necessary but not sufficient condition for this to occur is ii > 1 
or a baryon content more than 3 times the BBN value. 



4.5. Summary 

In summary, inflationary models carry an acoustic signature with distinct properties that can distinguish 
them from isocurvature models. The ratio of peak locations is a robust prediction of inflation. In particular, 
the ratio of third to first peak location should be in the range £3/^1 ~ 3.3 — 3.7 and the ratio of first 
peak location to peak spacing should be between Ii/Ia ~ 0.7 — 0.9. If this pattern is not observed in 
the CMB, either inflation does not provide the main source of perturbations in the early universe or BBN 
grossly underestimates the baryon fraction. The latter possibility is treated more fully in §5.6. Observational 
confirmation of the pattern would provide a plausibility proof for inflation. It would thus require fine tuning 
for an isocurvature model to reproduce this spectrum. To close this loophole in the test of inflation, the 
relative peak heights can be observed. Assuming the location of the peaks follows the inflationary prediction, 
we find that the high odd peak pattern of inflationary peaks is a unique prediction of inflation with a near 
BBN baryon content. Thus the locations and relative peak heights can be used to "prove" the inflationary 
paradigm. On the other hand, this signature relies on the baryon drag effect which may be small and 
difficult to observe in some exotic inflationary models. Moreover inflation does not preclude the presence of 
isocurvature perturbations (Linde 1985, Yokoyama & Suto 1991). Disproof of inflation along these lines is 
more difficult, a common problem in cosmology! 



5. Robustness of the Curvature Measurement 

Acoustic oscillations in the CMB provide an interesting and unique opportunity to measure the curvature 
of the universe. Features in the spectrum supply standard rulers with which to make a classical test of the 
geometry (Doroshkevich, Zel'dovich & Sunyaev 1978, Sugiyama & Gouda 1992, Kamionkowski, Spergel & 
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Fig. 11. — Standard rulers and the angular diameter distance. 

Acoustic features in the CMB, in particular the peak locations and the damping tail, act as standard rulers 

with which a measurement of the curvature can be made. Thin solid lines represent the Hq = 1, A = 
calculation scaled in I to account for the projection in the flo = 0.1 open and A models. 

Sugiyama 1994, Hu & White 1996). If the inflationary scenario is correct, the curvature can be measured to 
a few percent essentially from the location of the first acoustic peak. However, is the curvature measurement 
robust to changes in the underlying model? Clearly, the location of the first peak does not contain 
enough information. Isocurvature models predict peak locations with a different relation to the sound 
horizon. Radical changes in the baryon content from its BBN value and the thermal history from standard 
recombination can shift their positions. What additional information is required to ensure that the curvature 
measurement is robust? 



5.1. Angular Diameter Distance 

It is instructive first to review the general case for measuring the curvature by the angular diameter 

distance test. For dcfiniteness and simplicity, let us take the inflationary example. Here the acoustic spectrum 
in k is purely a function of VLbh? and Q-Qh"^ . The implied fluctuations on the last scattering surface are viewed 
as anisotropies on the sky today via the projection 

•^feature (-f'^) ■'^) = ^feature^S 

(K,A), (39) 

where the comoving angular diameter distance to the last scattering surface is 

re{K,K) = \K\-^'^ s\nh[\K\^'\'qo - rj^)]. (40) 

for ii' < negatively curved universes. For K > Q, merely replace sinh sin. The functional dependence 
arises since r]o{K,K,Q,Qh?) and ry*(f2o/i^, f^b/i^) and we are considering Q.oh? and Q.hh? fixed for now. In 
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general, negative curvature moves the acoustic features to smaller angles since a fixed physical scale then 
subtends a smaller angle on the sky. A cosmological constant moves features to larger angles due to a 
decrease in the current horizon size from the rapid expansion. 

In Fig. (11), we display an example. By holding f2(,/i^ and f2o/i^ constant, the acoustic features are fixed 
in physical space and the shift in I is entirely due to the projection. Shifting the canonical fio = 1-0 model by 
the angular diameter distance scaling of an fig = 0-1 model, we sec that Eq. (40) accurately accounts for the 
effect. Notice that K and A are degenerate with respect to the acoustic signature. One cannot in principle 
measure them independently from these features. The degeneracy is broken at larger angular scales by the 
ISW effect (HSb). Moreover for realistic Oq, A introduces little ambiguity in the curvature measurement. 

It is clear that either the peaks or the damping tail can be used as standard rulers to probe the curvature 
in this simple case. We shall now consider the benefits and drawbacks of each in the more general setting 
where the model for the perturbations and the background are not known a 'priori. We then turn to what 
can be learned about the model for the perturbations and how this additional information can be used to 
refine the measurement of the curvature. 



5.2. Acoustic Peaks 

The physical scale associated with the acoustic peaks is the sound horizon at last scattering. 
Unfortunately however, the exact relation between the peak locations and the sound horizon typically varies 
by a factor of 2 depending on the model. In this section, we discuss a how the peak spacing can be used as 
a measurement of the spatial curvature which is far less dependent on the model for the fluctuations. 

As discussed in §3.1, once the source has switched off inside the horizon the acoustic peaks follow a phase 
shifted harmonic series. The key point is that regardless of the phase, the peaks are spaced by kA = t^/ts- 
Hence if the sound horizon at last scattering is known, the spacing provides us with an angular diameter 
distance test of the curvature. In making this measurement, the higher peaks form a better probe than the 
lower peaks, since the latter can be shifted due to residual driving effects and metric fluctuations between 
last scattering and the present. The peak spacings A£ for the models of Fig. 7 are shown in Table 1 and 
should be compared with the prediction A£ = £a = kAVe ■ As can be seen, the spacing becomes more regular 
and approaches the expected value after the first peak separation ^2 — ^1- In the inflationary case shown here, 
baryon drag suppresses the even peaks and distorts their locations. For best acciiracy, one should employ 
the spacings of the prominent peaks for the test, e.g. £a = (^3 — ^i)/2 in the inflationary case. In general, 
three peaks will be necessary to assure accuracy of the test. 

How well can this test measure the curvature given realistic uncertainties in ilfe/i^ and Q.Qh? through h 
which affect the physical scale of the sound horizon at last scattering (see Fig. 12)? If the baryon content is 
near or less than the BBN limit Q.bh? < 0.05, it has only a small effect since the photons dominate the fluid 
at last scattering. The Hubble constant has a relatively larger effect but even so a measurement of a peak 
spacing £a ~ 290 would require f2o + ^ 0.7 if 0.4 <h< 0.8. Note that the dependence on the Hubble 
constant becomes weaker as it increases to make the universe more matter dominated at last scattering. 

The main drawback of this method is that it may be difficult to apply in models where the forced effects 
continue well after horizon crossing. In this case, the peak spacing may not become regular until the higher 
peaks. Near the diffusion scale, the power in the CMB fluctuations is exponentially damped. Foregrounds 
become more difficult to subtract on these smaller scales and gravitational lensing by large scale structure 
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(see e.g. Seljak 1996) can wash out the oscniations to some extent. It is even possible that stochastic metric 
perturbations would leave the spectrum with no distinct peaks (see §2.6, Magueijo et al. 1996). In these 
cases, we must turn to the damping tail to measure the curvature. 

5.3. Damping Tail 

The diffusion damping length provides a standard ruler for the curvature measiircmcnt that is the most 
robust against changes in the model for the fluctuations (sec also Hu & White 1996). As discussed in §3.4, 
its location is dependent only on background parameters. However, it is more sensitive to changes in r2(,/i^ 
and the thermal history due to its dependence on the Compton mean free path at last scattering. Fig. 12b 
displays the location of the damping tail £d defined as the multipolc number at which the acoustic effect 
drops by in power. It is related to the damping wavenumber Eq. (31) by the simple projection of 
Eq. (39). Here the thermal history has been fixed to follow standard recombination. By applying reasonable 
constraints on the other parameters, the damping scale becomes a sensitive probe of the curvature. Even 
allowing for a factor of 4 uncertainty in flbh"^, an open model of f2o ^ 0.5 can be distinguished from a flat 
fio + f^A = 1 universe. An additional uncertainty arises from the Hubble constant which changes the age 
of the universe and expansion rate at last scattering and hence the diffusion scale. In Fig. 12, we show the 
uncertainty from h if f2(,/i^ is fixed by BBN. Again open universes with Qq < 0.5 can be distinguished from 
the flat cases for reasonable values of h amounting to a factor of 4 ambiguity in floh^. 

There is one caveat to this proposal for measuring the curvature. It must first be established that 
the damping tail is indeed due to photon diffusion and not due to some intrinsic fall off in the source or 
secondary effect between recombination and the present. To test this assumption, one must measure the 
shape of the damping tail. Diffusion damping leads to a near exponential faU off in the anisotropy rather 
than the power law behavior expected in the alternate possibilities. Thus it will be necessary to follow the 
damped oscillations into the diffusion regime to establish its exponential character. Even in a flat universe 
the damping tail is at arcminute scales. Detailed measurements of the primary signal may or may not 
be observationally feasible due to possible foreground and secondary effects, most notably the ISW effect. 
Models with large small scale metric fluctuations may also have their effective temperature at last scattering 
dominated by baryon drag rather than acoustic oscillations. This too would mask the exponential signature. 
Of course, if acoustic peaks in the angular power spectrum are also measured, then it is almost assured that 
a break in the acoustic spectrum is due to diffusion. 



Model 


1-2 


2-3 


3-4 


4-5 


Inf 


288 


259 


297 


277 


HDM 


204 


282 


272 


295 


Tex 


221 


286 


287 


285 


Axi 


204 


295 


275 


298 



Table 1: The spacing between peaks for the suite of models in Fig. 7 with the same background parameters 
flo = 1, fife = 0.05, h = 0.5 and standard recombination. The spacing becomes more regular for the higher 
peaks and converges toward the expected value of 290. The spacing yields a measure of the curvature of the 
universe that is nearly independent of the model for the fluctuations. 
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Fig. 12. — Peak spacing and damping scale as a function of fio- 

Even allowing for uncertainties in the baryon content (solid shading, h = 0.5) and Hubble constant (dashed 
shading, Q^/i^ ^ 0.0125), open models with Qq ^ 0.5 can be distinguished from flat (fJo H-Ha = 1) A models 
through either scale. The damping scale is entirely independent of the model for the fluctuations but may 
be more difficult to measure than the peak spacing. 



5.4. Exotic Baryon Content 

Thus far, we have mainly considered the robustness of the curvature measurement to exotic sources 
of gravitational perturbations. It is also possible that the background itself exhibits an exotic nature. In 
particular, wc have been implicitly assuming that the baryon content is within a factor of several of its BBN 
value and the thermal history follows standard recombination. Let us now consider how the possibility of 
exotic background properties may be handled, beginning with the baryon content. 

If the baryon content is far from the BBN prediction of r2f,/i^ ~ 0.01-0.02, but known the approach to 
measuring the curvature is unchanged since both the sound horizon and damping scale are known as well. 
It may be the case however that the prediction is violated by observations that do not then constrain fl^h^ 
sufficiently. If it is low, the peak spacing still provides a good measurement of the curvature. In this limit, 
the sound horizon at a fixed redshift becomes independent of the baryon content and the only variation 
with Ofo/i^ is a weak dependence on the redshift of last scattering (see Fig. 13a). However, the peaks must 
still be measured for this test to work. With a low baryon content, the Compton mean free path and hence 
the diffusion length at last scattering increases making the higher peaks difficult to observe. The ratio of 
the damping tail to the peak-spacing Id/^-a provides a rough model independent estimate of the number of 
peaks that are potentially observable and is displayed in Fig. 13a (c.f. Fig. 5). Two peaks and hence the 
peak spacing is likely to be observable as long as n^h^ > 0.001, which is close to the stellar mass density 
f^stars ~ 0.004 for a reasonable Hubble constant. Thus the peak spacing should provide a test of the curvature 
even under the most extreme conditions. On the other hand, a curvature measurement from the damping 
tail alone coidcl fall victim to this exotic possibility. However, the severely truncated acoustic spectrum in 
this case should prevent such a misinterpretation of the position of the damping tail. 

In the case of an extremely high baryon content, the situation is reversed. The damping tail now provides 
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Fig. 13. — Exotic baryon and radiation content. 

(a) The sound horizon which forms the basis of the peak-spacing test is independent of O^/i^ for values near 
BBN or lower, whereas the damping tail is insensitive to high f2j,h^. In either limit a robust test exists. 
The ratio of tail to peak-spacing tojlA can detect an exotic baryon content, (b) Changing the radiation 
content, e.g. by altering the number of relativistic neutrinos Ni,, affects the expansion rate and thus the two 
physical scales weakly. Cases extreme enough to aifect the scales significantly can also be distinguished by 

the more robust estimate. Since the effects of a delay in last scattering and a decrease in the Compton mean 
free path tend to cancel, the physical scale is nearly independent of the baryon content for high values (see 
Fig. 13). On the other hand, if the baryon content is raised more than a factor of 5 over the BBN value, 
the sound speed and hence the sound horizon at last scattering decreases. Thus the peak spacing alone is 
an unreliable test of the curvature if the baryon content is extremely high > 0.06) but unknown. 

Is there a model-independent measure of the baryon content that can test for anomalous values? It 
should be clear from the discussion above that the ratio of the damping tail to peak-spacing is highly 

sensitive to Q,bh^ . It is also a fimction of fig^^ but since fig must be consistent with the inferred curvature 
only uncertainties in the Hubble constant enter if = 0. The ratio is double valued so that its measurement 
would allow both a high and low baryon content solution (see Fig. 13a). As we shall see in §5.6, baryon 
drag provides a distinction between the two extremes from a measurement of the relative heights of the 
peaks. In summary, an exotic baryon content is detectable and does not present a problem for the curvature 
measurement if both the peak spacing and damping tail can be measured. 

5.5. Exotic Thermal History 

Exotic thermal histories are another possibility. By thermal history, we refer to both the thermodynamics 

of the expansion and the ionization history. Massive decaying particles can create an epoch of matter 
domination before last scattering which changes the age of the universe at last scattering and hence both 
the sound horizon and the diffusion scale (White, Gelmini & Silk 1995, Bardeen, Bond & Efstathiou 1987). 
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Fig. 14. — Exotic ionization history 

We show the damping and peak-separation scales in a model with instantaneous recombination at z* . For a 

gradual recombination, a lower lol^-A will always result. The diffusion scale at last scattering by definition 
approaches the sound horizon in the limit that no recombination (NR) occurred. Assuming instantaneous 
recombination, the ratio lo/tA roughly corresponds to the number of observable peaks and can be used to 
discriminate against exotic ionization histories. Here Vloh? = 0.25, il^h^ = 0.0125. 

A less exotic example of the latter effect is provided by any model which changes the epoch of equality. 

Let us take the a simple but illustrative example of a change in the energy density in relativistic species. 
This can arise for example by changing the number of relativistic neutrinos or their temperature. In Fig. 13b, 
we plot the sound horizon, damping scale and the ratio between the two as a function of the fractional increase 
in the radiation density Pr/ p-y- The standard thermal history with three relativistic neutrinos predicts 
Pr/ P-y = 1-68. Raising the radiation content increases the expansion rate and thus decreases both the sound 
horizon and diffusion scale at last scattering. Its effects are relatively weak unless the standard prediction 
grossly underestimates the radiation content. Because the diffusion scale is essentially the geometric mean 
of the horizon and Compton mean free path, it is a weaker function of the radiation content. Hence, if the 
universe really has a sufficiently exotic history to change the two physical scales at last scattering, it should 
be detectable in the ratio of tail location to peak spacing loji-A- 

The use oi luj^A as a probe of the thermal history is especially powerful for models which delay last 
scattering significantly. This might occur due to energy injection near z ^ 10"^ from particle decays or 
non-linear fluctuations. Since the Compton mean free path grows as the universe expands, the diffusion 
scale can approach the horizon scale (see Fig. 14). Here full ionization is assumed until a redshift z^,, when 
recombination is taken to occur instantaneously. Instant recombination is not realistic; however it provides 
an upper limit on £o- The ratio of £d/C-a in Fig. 14 then measures the delay in recombination. Under 
the instantaneous recombination assumption it corresponds to the number of observable peaks. The actual 
observable value for any given model will be lowered in a manner dependent on the details of recombination. 



-35- 



Even if recombination never occured the Compton mean free path will eventually reach the horizon size 
due to dilution in the electron number density from the expansion. At this point the universe is effectively 
transparent and the photons free stream to the present. Since by definition the diffusion length roughly 
corresponds to the horizon scale at this epoch, IdI^a ~ 1 and no acoustic oscillations are apparent. Thus 
recombination is a necessary condition for the acoustic signature to be observable. In fact for practical 
purposes, we require at least two observable peaks. For fio/i^ = 0.25 and = 0.0125, this translates 

to recombination by = 175. For comparison, late reionization after standard recombination does not 
destroy the acoustic signal unless the optical depth through the reionized epoch is greater than unity: 

^reion Z 100(Oo/lV0.25) V3(J]j/i2/q_0125)-2/3. 

Although there is not enough information in the one number I-q / (.a to reconstruct all possible exotic 
thermal histories, it is possible to determine that some exotic model is necessary if IbI^a is not « 4 — 5. 
Further work on the damping tail is necessary before this constraint can be tightened (Hu & White, in 
preparation). In this case, we trade precise knowledge of the curvature for evidence that exotic physics is 
required in the early universe. On the other hand, if the deviation is known to occur due to a specific cause 
such as early energy injection or delays in equality, then the curvature can once again be obtained. Moreover, 
with this information from both the damping tail and the peaks, we cannot mistakenly infer a value of the 
curvature due to exotic baryon content or thermal history. 

5.6. Toward Reconstructing the Model 

The previous sections considered constraints on the curvature that could be made with minimal 
knowledge of the model for the perturbations. However we would ideally like to learn as much as possible 
about both questions from the CMB spectrum. In this section, we shall outline a program for measuring the 

curvature by first establishing the basic model for structure formation. Model properties such as the relative 
peak heights simplify the task of measuring the curvature since they provide extra clues. As noted in §§5.2 
- 5.3, this may be the only way to proceed if it proves impossible to measure the higher acoustic peaks and 
damping scale. 

We define the basic model possibilities by their dominant mechanism for forming anisotropics : 

1. Inflation: correlated superhorizon curvature fluctuations. Although inflation may also produce 
isocurvature initial conditions, we class them separately. 

2. Driven isocurvature (e.g. axionic and baryonic isocurvature, textures, etc.): photon-compensated 
initial conditions. We allow for the possibility that the first peak is obscured (see §4.2) especially in 
high baryon cases where it is strongly suppressed with respect to the second peak. 

3. Forced isocurvature (e.g. possibly strings): forces after but temporally correlated with horizon 
crossing. 

4. Stochastic isocurvature (e.g. possibly strings): random sub-horizon forces. 

5. Reionized: any of the above in which the diffusion length and horizon length coincide, e.g. fully 
ionized models. 

Each of these models can have anomalously high (low) baryon content defined as Ofo/i^ significantly larger 
(smaller) than the big bang nucleosynthesis (BBN) value of 1-2%. They can also be low in dark matter 
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(CDM) such that fig ~ ^if,. It is possible, though unUkely, that the model may be a mixture of "inflation" 
(as defined above) and driven isocurvature scenarios. 

For dcfinitcncss, wc hereafter discount the possibility of an exotic thermal history as treated in §5.5, 
except for the case of late reionization which is plausible in many models. Furthermore, we assume that the 
Hubble constant and A are suflaciently constrained to make the confusion they introduce to the curvature 
measurement irrelevant (see e.g. Fig. 12). We however relax the BBN constraints on f2(,/i^ as this can lead 
to a qualitatively distinct acoustic pattern. 

There are four tests that we can apply to fix the model and the curvature, based on the acoustic 
signatures discussed in §3. 

1. Peak positions test: measure the locations and spacings between the peaks in I. Three useful items 
can be extracted from this test: 

Peak ratios: Ratios probe the nature of the model. We define the distinguishing feature of a "cosine" 
series as a ratio of third to first peaks of I'^jlx « 3 — 4 and that of a "sine" of « 5. Note that a true 
sine oscillation with missing first peak might be observationally classed as a cosine. 

Peak spacings: If they are regular, they provide an angular diameter distance test. The physical scale 
corresponds to kA = Tr/r^ if both compressional and rarefactional peaks are measured or 27r/rs for the 
compressional ones. 

Peak-to-spacing ratio: Since the peak spacing is fixed by the sound horizon, its ratio with respect to 
the first peak location provides a sensitive probe of the model. Inflation predicts ixjiA ~ 0.7 — 0.9. 

2. Relative height test: determine if any of the peaks are anomalously high with respect to a smooth 
underlying spectrum. Heights probe the baryon content and forcing mechanism. 

3. Tail-to-spacing test: measure the ratio of tail to peak-spacing Id/^a- If the peak spacing is regular, 
the ratio measures the baryon content Q-bh? (and/or identifies exotic thermal histories). 

4. Damping tail test: measure the shape and absolute location of the damping tail as described in §5.3. 
The shape confirms its acoustic nature and the location provides an angular diameter distance test. 

As we shall sec, in most cases not all four tests are necessary. In contrast to §5.3, we shall here adopt the 
philosophy that information from the damping tail only be used if all other tests are ambiguous. If the 
baryon content is assumed to be known from nucleosynthesis, most of the tests involving the damping tail 
are generally unnecessary. On the other hand, the damping scale ideally should still be measured since it 
provides a valuable consistency check on both the baryon content and the thermal history assumption. Let us 
outline the program as a fiow chart. It starts with the suite of peak position tests. There are 4 possibilities: 

A. The peak ratios form a cosine series. 

Apply the relative heights test. 

1. High odd peaks forming a smooth sequence: inflation with BBN or high Q.^h?'. Tail-to-spacing fixes 
fifo/i^. Spacing measures the curvature. 

2. High odd peaks with anomalously high flrst peak: high Q,i,h? driven isocurvature with first peak 
obscured. Tail- to-spacing fixes O^/i^. Spacing measures curvature. 
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3. High even peaks: tuned forced isocurvature with near BBN Spacing measures the curvature 

4. Monotonic: BBN or lower Qbh"^. Spacing measures curvature. 

a. Peak-to-spacing ratio < 1: low CDM inflation. Tail- to-spacing can also measure ilih^. 

b. Peak-to-spacing ratio > 1: driven isocurvature with obscured first peak. Tail-to-spacing can also 
measure O^/i^. 

B. The peak ratios form a sine series. 

Apply the relative heights test. 

1. High even peaks: driven isocurvature model with high Cl^h^. Tail-to-spacing fixes flhh^. Spacing or 
damping tail measures curvature. 

2. High second peak: driven isocurvature model with BBN or lower. Spacing measures curvature. 
Tail-to-spacing can also measure flhh^. 

3. Monotonic: high fl^h^, high CDM inflation. Tail-to-spacing fixes Oft/i^ and /la/f spacing determines 
curvature. 

C. Neither cosine nor sine 

1. Peaks follow phase shifted harmonic: Coherent and tuned superposition of inflationary and driven 
isocurvature (cf. §3.1). Tail-to-spacing fixes flhh^. Spacing measures curvature. 

2. Peaks follow a cosine harmonic with gaps 1:3:4:5: High Qbh"^, relatively high ^o/^b inflation. 
Tail-to-spacing fixes flbh^. Spacing of higher peaks or /ia/!fspacing across the gap measures curvature. 
Damping tail also measures curvature. 

3. First few peaks irregular followed by a regular series. Forced isocurvature model with horizon crossing 
effects and BBN or lower f2(,/i^. Peak spacing measures curvature. Tail-to-spacing ratio can also 
measure the 

4. Somewhat irregular peaks offset by smooth function (e.g. Fig. 3). Forced isocurvature model 
with effects well after horizon crossing and high flbh^. Take average value of half-peak spacing. 
Tail-to- (average)-spacing fixes ftbh^. Average spacing measures curvature. 

5. Random locations: Forced isocurvature model with rapidly varying metric well after horizon crossing. 
Tail-to-spacing ratio constrains Ofe/i^. Damping tail constrains curvature. 

D. No peaks. 

Damping tail test. 

1. Exponential fall off: Stochastic isocurvature model. Location of damping tail measures the curvature 
if Clbh^ is known. 

2. Power law fall off: Reionized or stochastic isocurvature model with peak source amplitude far inside 
horizon and high No robust constraints on the curvature are possible. 
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Thus, in all but the last case the acoustic signature constrains the curvature. If the baryon content can 
also either be measured from the signature itself or is known from external constraints such as BBN, highly 
accurate measurements of the curvature are possible. Once the basic nature of the model and background 
is determined through this program, detailed modeling of the fundamental source for the fluctuations that 
formed large scale structure may begin. 



6. Conclusions 

We have generalized the formalism of HSa to include backreaction effects and examined the uniqueness 
and robustness of acoustic signatures in the CMB. By clarifying the role of compensation and feedback in 
the evolution of fluctuations, we have shown that the phase of the oscillation, and hence the ratios of peak 
locations, distinguishes inflation from "typical" isocurvature models. Specifically, two robust test are the 
ratio of third to first peaks and first to peak-spacing. Our analysis also provides a better understanding of 
the structure of the peaks in these models. However since it is possible to imagine isocurvature models which 
are tuned to mimic the inflationary pattern of peaks, we have stressed the importance of baryon drag, which 
allows us to distinguish compressions from rarefactions in potential wells. This can help lift the confusion 
between adiabatic models and these contrived isocurvature models. Although the level of drag does not 
make this distinction clear for all possible baryon densities, it is observable for the value predicted by big 
bang nucleosynthesis. Further understanding of the diffusion damping of anisotropics will allow us to extend 
the lower limit of baryon densities for which we can distinguish the rarefaction and compression peaks. 

We focus on the importance of the damping tail as a measure of spatial curvature which is independent 
of the model for structure formation, and discuss the robustness of curvature measurements from the CMB. 
Even in the case where no acoustic peaks are seen (e.g. possibly string models) the damping scale can be 
estimated under mild assumptions about the thermal history and baryon content. A more precise test is 
possible if the acoustic peaks are regularly spaced as indeed expected of models without extreme behavior at 
small scales. Either approach allows one to infer the physical scale of the acoustic feature(s). Its projection 
onto the sky allows us to perform a classical angular diameter distance test to determine the curvature of 
the universe. In the process of carrying out these tests of the curvature, the general nature of the model for 
the fluctuations can be reconstructed as well as the baryon content of the universe. 

All of these studies focus on the tale told by the CMB spectrum taken as a whole. In particular, the 
acoustic pattern, which arise from forced oscillations in the photon-baryon fluid before recombination, leaves 
a distinct signature from which we may begin to reconstruct the cosmological mdoel. 
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A. Causality and Compensation 

In this Appendix, we clarify the role of causality in limiting the behavior of fluctuations outside the 
horizon and its dependence on gauge. We show that only inflation can correlate the curvature fluctuations 
above the horizon assuming general relativity is the correct description of gravity. Moreover in the main 
text, we employed only the mechanism of compensation, i.e. the response of the photon-baryon fluid to 
a source, and not the full causal constraints which limit the behavior of the source as well. This proved 
sufficiently powerful to produce the distinctions in the acoustic signature under the additional assumption 
that the compensation is provided by the photons. In this case, the feedback from the photon self-gravity 
produces the key ingredient in making many of the signatures robust. Since in the standard FRW model, 
the universe is radiation dominated until near recombination, this additional assumption is automatically 
satisfied. There arc however exotic models where this assumption is not satisfied. For example, a decaying 
massive particle could cause the universe to undergo a period of matter domination before recombination. For 
this kind of situation, we need to examine the general case of compensation and additional causal constraints 
on the model. 

Causality implies initial compensation in density fluctuations above the horizon since the stress-energy 
tensor is conserved. Heuristically the conservation law implies that changes in energy density at any location 
arise from "flows" of energy density current across surfaces, or from displacements of fluid elements. Since 
fluid elements cannot be displaced "beyond" the horizon, this severely constrains the behavior of fluctuations 
at fcr/ <C 1. However, because the stress-energy tensor obeys covariant conservation, the exact form of causal 
constraints on the density perturbation depends on the representation of the metric, i.e. the gauge. 

We have seen in Eq. (1) that density perturbations can also change due to the "stretching" effects from 
changes in the metric. A clever choice of gauge can eliminate such effects. Let us examine the evolution of 
the density fluctuation in its local rest frame. Note that this does not coincide with hypersurfaces with zero 
bulk velocity for the total matter Vr unless goi vanishes. In the literature, this has been called the comoving 
gauge (Bardecn 1980), velocity-orthogonal isotropic gauge (Kodama & Sasaki 1984) and total matter gauge 
(HSb). It represents the metric fluctuations as 

500 = -aHl + 2^Q), 

goj = a''VTk-'Q\j, (Al) 
Qij = a^{l + 2CQ)jij, 

where | represents a covariant derivative with respect to ^ij and recall Q is the kth eigenfunction of the 
normal mode decomposition. These quantities are related to their Newtonian counterparts by a gauge 
transformation x^^ = + Sx^^. If the line element ds"^ is to remain unchanged (to first order in Sx), 

Qnv = 9tiu + gai^Sx" ^ + gaf_iSx" ,, - gf,„,aSx". (A2) 

A similar relation follows for the stress energy tensor T^^ and relates the matter quantities. For densities, 
it is simpler to note that it arises from a combination of a shift in time slicing and the background density 
evolution: SiQ = SiQ — {pi/pi)Sx°. These relations imply that from the Newtonian gauge, the rest frame is 



'^Note that the universe need only be radiation dominated when the fluctuation was well outside the horizon, not at last 
scattering itself. This is of course not satisfied in models with significant reionization, but in this case ax;oustic oscillations are 
not observable anyway. 
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reached by the coordinate shift 6x^ = {Vr/k, 0)Q, and the perturbation quantities are related by 

<^ = $ - -Vr/k, (A3) 
a 

At = ST + 3{l + WT)-VT/k, 
a 

where Wi = Pi/pi and subscript T denotes quantities of the total fluid from a sum over the particle 

constituents. Here wc have used the background evolution equation pi/ pi = —3(1 + Wi){a/a). Note that the 
velocity Vr is the same in the Newtonian and rest frame gauges since its transformation properties depend 
on 6x'^. The familiar quantity ^ is the curvature of the spatial hypersurfaces in this gauge (Bardeen 1980). 

Applying the gauge transformation Eq. (A3) to the Newtonian gauge equations or writing the Einstein 
equations in this gauge (see HSb, Eq. (16), Kodama & Sasaki 1984 Eq. 4.7), we obtain the evolution equation 
for the curvature, 

^ al + WT \Pt S J ^ ' 

Unless wt = —1 as in the case of the de Sitter phase, the rest frame curvature remains constant in the 
absence of isotropic (pressure) or anisotropic stress perturbations. The isotropic stress perturbation can be 
broken up into an adiabatic and non-adiabatic part 

^ =(^w^^At + Tt, (A5) 
Pt 

with the sound speed of the total fluid given by = Pt/ Pt- The adiabatic pressure perturbation is related 
to the curvature fluctuation by a factor of (fcry)^ through the Poisson equation and hence is a negligible source 
for [kf])'^ <C 1. The generalized compensation law is therefore that component evolution must balance at 

krj to keep the rest frame curvature constant in the absence of non-adiabatic isotropic stress ( "entropy" ) 
or anisotropic stress perturbations. Note that unlike the adiabatic isotropic stress, both of these sources are 
the same in any frame and admit no gauge ambiguity. 

Since the continuity equation in this gauge is 

for each of the individual particle species, the rest frame density perturbation obeys an ordinary conservation 
law if these stresses are absent (C = 0). In other words, the number density of each of the particle constituents 
only changes through their bulk motion. Thus we see that the causal constraint is simplest in the rest frame. 
Note that Eqs. (A6) and (A4) imply that although density fluctuations can change purely due to an evolving 
equation of state, this does not affect the curvature fluctuation C, above the horizon. 

Now let us consider the source of curvature fluctuations from non-adiabatic pressure perturbations. The 
analysis also applies to anisotropic stress perturbations. The curvature fluctuation generated by pressure 

perturbations is of order 5pT / (pt + Pt) ■ This fact is somewhat counterintuitive since the physical mechanism 
that converts a pressure fluctuation to a density fluctuation is the movement of matter which is impossible 
beyond the horizon. Let us examine its qualitative origin. Causality constraints in Fourier space do not 
require no evolution for kr] <^ 1. In physical space, motion of matter through length scales up to rj cause a 
suppressed evolution of the Fourier amplitude. A change in momentum density of the fluid is caused by the 



(A6) 
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pressure gradient and generates a bulk velocity of order [kifjSpT / {pt + Pt)- This then forms a kinematic 
density fluctuation from the continuity equation (A6) of order {krj)'^ Spx / {pt + Pt) or a curvature fluctuation 
of order Spx/ipr + Pt)- Thus the residual curvature fluctuation induced by motion of matter inside the 
horizon is generically of order Wt/(1 + Wt) times the pressure perturbation. 

Causality implies that before matter has had a chance to move around, the universe must obey the 
isocurvature condition ^ = or At =0.^ If the condition is met by balancing perturbations in the different 
components of the fluid, a non-adiabatic pressure fluctuation generically arises 

PtTt = SpT - c%8pT = ^(c^ - 4)(5pi, (A7) 

if the equation of state of the balancing components differ. As we have seen, can produce a curvature 
fluctuation even at <C 1. How docs this possibility affect our arguments concerning the uniqueness 
of the inflationary spectrum? First, we need to relate the rest frame curvature to the Newtonian 
curvature. Employing the continuity equation (Kodama & Sasaki 1984, Eq. 4.7) in Eq. (A3), we find 
(Lyth 1985, Mukhanov, Feldman & Brandenberger 1992) 

C = * + |t^(* + t^), (A8) 
31 + wt ^ a J 

in the absence of anisotropic stress. If the equation of state is constant and $ evolves as a power law, C oc 
The two curvature fluctuations are comparable except in the degenerate case where C = and 

| = -^[^(1 + «'t) + 1]. (A9) 

In fact, this is the case of the source-free decaying mode described by Eq. (14). The decaying mode thus 
carries no curvature perturbation in the rest frame. The other case in which and $ behave differently is 
through a change in the equation of state. For example, through the matter radiation transition wt goes 
from 1/3 0. Although ( remains constant in the absence of stress perturbations, $ drops by a factor of 
9/10 through the transition. For most purposes however, we can think of and $ as nearly interchangeable. 

The mechanism by which non-adiabatic pressure perturbations generate curvature fluctuations is of 
course already implicitly encorporated in our analysis and is the cause of the curvature not being strictly 
zero outside at kr] <^ 1 but merely suppressed. In the photon dominated limit, only a small pressure 
fluctuation is needed to compensate a rather large density fluctuation in the source. Thus, the curvature 
generated by this effect is negligible. In the matter dominated limit, pressure fluctuations cannot move a 
large amount of energy density as exhibited by the wt/{1 + wt) suppression. However, as the universe 
changes from radiation to matter domination a relatively significant curvature fluctuation can be generated. 

The baryon isocurvature model provides a concrete example. In this case, the non-adiabatic pressure 
perturbation is 

rT = -i^5, (AlO) 

where recall that S = 6{ni,/n^) = Sb — 39o is the entropy fluctuation. Notice that the pressure fluctuation 

is small as long as the universe is radiation dominated wt = 1/3. The evolution equations may be exactly 
solved in the krj limit such that the Newtonian curvature is $ = ^{a/aeq)S in the radiation-dominated 
limit and in the matter dominated limit (see HSb, Eq. 27). Notice that the curvature is a constant in 



^More precisely Aj- = aside from a C = decaying mode [see Eq. (A9)]. 
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the matter dominated limit even outside the horizon. Around equality, pressure perturbations of order S 
generate curvature fluctuations of the same order. Since in the matter dominated limit, pressure is no longer 
effective, this curvature fluctuation is frozen in. 

If one considers the evolution of a single /c-mode, distinguishing between this and the inflationary 
case would be difficult since they both exhibit a constant curvature fluctuation above the horizon. In the 
standard scenario, this is not a problem since the matter-radiation transition cannot occur early enough 
without overdosing the universe (ilo^^ ^ !)• Note that the curvature fluctuation must be constant well 
before horizon crossing for all observable scales in order to mimic inflation. However it is possible with 
decaying massive particle scenarios for the universe to undergo a period of matter- domination before the 
ordinary radiation-dominated epoch. 

In this case, the spectrum of non-adiabatic pressure perturbations implied by causality serves to 
distinguish the model from inflation. Causality forbids spatial correlations in the spectrum of such 
perturbations, so |rT(fc)p and hence are constant in k, i.e. white noise (or steeper if symmetries 

can be imposed), before horizon crossing. Notice that this agrees with the familiar result that the density 
perturbations generated by the motion of matter have a | A^p cx tail for kr] <C f (Zcl'dovich 1965, Robinson 
& Wandelt 1996). This translates to a steeply rising spectrum of acoustic fluctuations compared with the 
inflationary case of an approximately scale invariant spectrum |$(?7j,/c)p oc k~^. Thus acoustic modes 
associated with a constant curvature perturbation outside the horizon generated from an isocurvature initial 
condition are both easily distinguished from inflation and observationally ruled out! The most general 
isocurvature spectrum with scale invariant curvature perturbations at horizon crossing is |<I'(?y, fc)p = 
F{kri)k~^. With the requirement of white noise perturbations outside the horizon, $ must grow as rf^'^ 
before horizon crossing. Since only models like these, with curvature fluctuations growing until horizon 
crossing, require observation of the acoustic signature to distinguish them from inflation, our assumption in 
the main text is justifled. It is of course still possible that tuned effects around horizon crossing can mimic 
an inflationary spectrum in such a exotic scenario where the regulatory effects of photon feedback are absent. 
However, since both tuning and a drastic modification of the thermal history is necessary, we do not consider 
this possibility to be worth considering. 
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